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ABSTRACT 


An  iterative  numerical  technique  for  analysing  the  biaxial  res¬ 
ponse  of  reinforced  and  prestressed  concrete  wall  segments  subject  to 
combinations  of  prestressing,  creep,  temperature  and  live  loads  is  pre¬ 
sented  . 

Two  concrete  constitutive  relations  are  available  for  this  analysis. 
The  first  is  a  uniaxially  bilinear  model  with  a  tension  cut-off.  The 
second  is  a  nonlinear  biaxial  relation  incorporating  equivalent  uniaxial 
strains  to  remove  the  Poisson’s  ratio  effect  under  biaxial  loading. 

Predictions  from  both  the  bilinear  and  nonlinear  model  are  compared 
with  observations  from  experimental  wall  segments  tested  in  tension. 

The  nonlinear  model  results  are  shown  to  be  close  to  those  of  the  test 
segments,  while  the  bilinear  results  are  good  up  to  cracking. 

Further  comparisons  are  made  between  the  nonlinear  analysis  using 
constant  membrane  force-moment  ratios,  constant  membrane  force-curvature 
ratios,  and  a  nonlinear  finite  difference  analysis  of  a  test  containment 
structure.  Neither  nonlinear  analysis  could  predict  the  response  of 
every  wall  segment  within  the  structure,  but  the  constant  membrane  force- 
moment  analysis  provided  lower  bound  results. 
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1.  INTRODUCTION 


1.1  Background 

During  the  last  five  years,  the  Department  of  Civil  Engineering 
at  the  University  of  Alberta  has  been  conducting  a  research  program 
to  aid  in  the  understanding  of  the  overload  behavior  and  design  pro¬ 
blems  associated  with  nuclear  secondary  containment  structures.  The 
structure  is  modelled  as  a  segmented  axisymmetric  shell  constructed  of 
concrete,  reinforced  and  prestressed  in  two  directions.  Due  to  the 
prestressing  and  the  load  carrying  mechanisms  of  a  shell,  the  walls 
and  roof  of  the  containment  structure  are  subject  to  in-plane  normal 
forces  and  out-of-plane  moments  in  both  directions  simultaneously. 

Under  the  overloads  to  which  the  walls  and  dome  of  the  structure 
may  be  subjected,  a  large  portion  of  the  vessel  would  be  stressed  in 
tension,  uniaxial  or  biaxial.  This  deviates  from  the  usual  stress 
state  to  which  concrete  is  subjected  at  design  loads  and,  therefore, 
the  familiar  assumptions  and  theory  applied  to  concrete  in  tension  may 
not  give  a  realistic  assessment  of  the  behavior  of  the  containment. 

A  previous  study  addressing  the  problem  of  analysing  a  wall  seg¬ 
ment  from  a  containment  vessel  was  conducted  by  Epstein  and  Murray 
(1976).  In  addition  to  computing  the  forces  acting  at  design  load  on 
the  wall  segments  variously  located  through  the  structure,  Epstein  and 
Murray  (1976)  attempted  to  establish  load  factors  to  various  limit  states 
by  analysing  the  behavior  of  these  segments.  This  analysis  was  elastic 
and  one-dimensional,  and  considered  the  transformed  area  and  moment  of 
inertia,  with  due  regard  to  cracking,  to  calculate  the  wall  segment 
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response.  The  behavior  of  concrete  in  tension  was  assumed  to  follow 
the  classical  theory  of  being  linearly  elastic  up  to  its  fracture 
strength  and  to  lose  all  stiffness  after  cracking.  (This  will  be 
referred  to  as  a  "tension  cut-off"  analysis) . 

1.2  Objectives  and  Scope 

The  shell  element  to  be  analysed  herein  is  similar  to  that  dis¬ 
cussed  in  the  report  of  Epstein  and  Murray  (1976)  and  is  depicted  in 
Fig.  1.1.  The  element  is  of  arbitrary  size,  doubly  curved  and  has  two 
degrees  of  freedom  in  each  of  the  two  orthogonal  directions  associated 
with  the  principal  radii  of  curvature.  These  degrees  of  freedom  allow 
in-plane  stretching  deformations  and  out-of-plane  curvatures.  The 
analysis  developed  herein  will  accept  the  presence  of  reinforcing  steel 
or  prestressing  steel  running  in  the  direction  of  either  of  the  princi¬ 
pal  curvatures.  These  steel  elements  can  be  placed  at  any  location 
within  the  thickness  dimension  of  the  wall. 

The  objective  is  to  develop  a  computer-aided  numerical  analysis 
to  predict  the  biaxial  response  of  wall  segments,  using  either  a  uni¬ 
axial  bilinear  elastic  concrete  constitutive  relation,  or  a  biaxial 
nonlinear  concrete  constitutive  relation,  and  to  compare  the  results 
obtained  from  this  program,  which  shall  be  referred  to  by  the  acronym 
COWSAC  (COWSAC  is  a  short  form  of  Containment  Wall  S_egment  Analysis  Code) , 
with  those  from  experimental  tests  of  wall  segments  and  from  other  non¬ 
linear  analyses. 

Two  concrete  constitutive  relations  will  be  incorporated  into  the 
analysis.  The  first  will  be  a  bilinear  elastic  model  with  a  tension 
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cut-off  considering  each  dimension  independently,  similar  to  that 
employed  by  Epstein  and  Murray  (1976).  The  second  will  be  a  two-dimen¬ 
sional  nonlinear  analysis  including  the  softening  behavior  of  the  con¬ 
crete  stress-strain  curve  in  both  tension  and  compression.  This 
unorthodox  treatment  of  the  behavior  of  concrete  in  tension  is  expected 
to  allow  simulation  of  problems  where  the  wall  segment  is  loaded  in 
tension.  The  prestressing  and  reinforcing  steel  will  be  modelled  with 
uniaxial  piecewise  linear  stress-strain  relations. 

In  the  present  work,  the  "load"  applied  to  the  wall  segment 
associated  with  each  degree  of  freedom  can  be  either  a  force  or  a  defor¬ 
mation.  Thus,  for  the  axial  degrees  of  freedom,  either  a  normal  force, 

Nj  or  N2,  or  a  normal  strain  or  e2  >  can  be  imposed.  For  the  rota¬ 
tional  degrees  of  freedom,  either  a  moment,  or  M2,  or  a  change  in 

curvature,  <j>]_  or  <j>2,  may  be  imposed.  With  a  total  of  four  degrees  of 
freedom,  and  a  choice  of  two  types  of  "loads"  associated  with  each, 
sixteen  possible  types  of  loading  patterns  may  be  specified. 

In  addition  to  the  live  loads,  the  wall  segment  may  be  subject  to 
prestressing  loads,  creep  effects,  and  thermal  loads.  For  each  of  these 
load  states,  the  analysis  can  account,  through  the  strain  displacement 
equations,  for  nonuniform  strains  across  the  section. 
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Fig.  1.1 


Wall  Segment  Showing  Degrees  of  Freedom 


2 .  THEORY 


2.1  Introduction 

The  fundamental  theory  used  to  model  the  wall  segment  is  presented 
in  this  chapter.  The  assumptions  which  pertain  to  both  the  material 
stress-strain  relationships  of  the  steel  and  concrete,  and  the  wall  seg¬ 
ment  theory  are  first  presented.  The  general  theory  for  a  wall  segment 
is  then  discussed.  Central  to  this  theory  are  the  wall  segment  equilib¬ 
rium  equations  developed  to  include  the  segment  layering,  the  interaction 
of  the  layers,  and  the  response  of  the  segment  to  either  imposed  forces 
or  strains.  The  Newton-Raphson  solution  technique  is  then  discussed 
relative  to  these  equations.  Finally,  the  specialized  problems  of 
thermal  and  creep  loads  are  treated. 

2.2  Assumptions 

To  simplify  the  mathematical  formulation  of  the  problem,  the 
following  assumptions  regarding  the  plane-stress  analysis  of  the  wall 
segment  are  made. 

1.  The  Love-Kirchhof f  assumptions  are  valid.  That  is,  normals 
to  the  reference  surface  remain  straight  and  normal  after 
deformation.  Material  layers  can  be  considered  to  be  in 

a  state  of  plane  stress. 

2.  No  in-plane  shear  forces  are  present.  That  is,  principal 
stresses  and  principal  directions  always  coincide  with  the 
directions  of  principal  curvature.  This  condition  occurs 
when  axisymmetric  structures  are  loaded  axisymmetrically . 

3.  Strain  comparability  is  always  maintained,  in  both  elastic 
and  inelastic  regions,  between  the  concrete  and  steel. 
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4.  The  steel  layers  carry  uniaxial  stress  only. 

5.  The  steel  layers  may  be  modelled  as  continuous  layers  of 
equivalent  area  extending  across  the  width  of  the  segment. 
The  size  and  spacing  of  the  steel  bars  or  tendons  have  no 
effect  other  than  determining  the  effective  thickness  of 
the  layer. 

2.3  Fundamental  Theory 

The  development  of  equations  parallels  that  for  B0S0R5  (Bushnell, 
1974)  and  is  based  on  the  principle  of  virtual  work.  The  program 
simulates  the  relation  between  the  force  vector  {Q}  and  the  correspond¬ 
ing  displacement  vector  {q},  where 

{Q}  =  <N]_ ,  Mi,  N2,  M2>  (2.3.1) 

and 

{q}  =  <£imS,  4> i  ,  t2mS>  <t>2>  (2.3.2) 

These  forces  and  deformations  are  shown  in  Fig.  1.1. 

The  object  of  the  program  is,  given  a  force,  to  find  the  resulting 
strain  or,  given  a  strain,  to  find  the  corresponding  force,  so  that  at 
the  time  of  solution,  both  {Q}  and  {q}  are  known. 

From  the  principle  of  virtual  work 

6U  =  a.  6e  dV  =  Q .  6q .  =  6W  (2.3.3) 

V  k  k  J  J 

where  6U  represents  the  virtual  strain  energy,  V  is  the  volume  of  the 

body,  a,  represents  a  vector  (consisting  of  the  stresses  0\  and  a2), 

e  represents  a  vector  consisting  of  the  corresponding  strains  ej  and 
lc 

e2,  6W  is  the  virtual  work  done  by  the  external  forces  and  6  is  the 
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standard  symbol  denoting  a  change  in  the  variable  which  it  precedes. 
Thus  the  summation  convention  is  used  in  Eq .  2.3.3,  with  a  range  of  2 
for  k,  and  a  range  of  4  for  j. 

If  the  matrix,  D.  .  is  now  introduced  which  relates  a,  and  e.  by 

ki  k  l 

the  following  equation 


then  Eq. 


2.3.3 


£  . 
1 


can  be  written  as 


<5U 


V  Dki 


£ .  5e.  dV 

l  k 


Q  .  6  q  . 
J  3 


6W 


(2.3.4) 


(2.3.5) 


Equilibrium  requires  that  the  external  forces  on  a  body  must  be 
resisted  by  equal  forces  developed  within  the  body.  By  defining  T ,  a 
variable  force  which  represents  the  difference  between  the  internal  and 
external  forces,  one  can  write 


r 


i\^dV-Q  -  0 


(2.3.6) 


Since  the  internal  and  external  forces  are  equal,  ¥  should  always  equal 
zero  in  an  equilibrium  condition. 


An  increment  in  T  can  be  represented  by 


AT, 


3T.  .  n 

_J.  K-  =  0 


8qi 


or  from  Eq.  2.3.6  by 


(2.3.7) 


AT.  = 
J 


/ 

V 


9cr,  9e, 


2 

9  e. 


9q .  9q  . 
J 


+ 


9q  .  9q  . 
i  J 


dVAq^  -  AQ_. 


=  0  (2.3.8) 


' 
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If  the  unbalanced  forces  Y.  are  not  zero  then  corrections  to  q. 

J  i 

can  be  obtained  by  requiring 


AY.  =  -  Y . 

J  J 


(2.3.9) 


Equations  2.3.9  are  simultaneous  non-linear  equations  which  are  known 
as  the  Newton-Raphson  equations.  Using  Eq .  2.3.7,  Eq .  2.3.9  may  be 


rewritten  as 


.1  Aq  =  -Y  (2.3.10) 

3t>i  J 

where  the  subscripts  i  and  j  have  a  range  of  N,  where  N  equals  the  number 
of  degrees  of  freedom,  which  in  this  case  is  four. 

The  above  equations  are  applied  to  a  wall  segment  which  is  layered 
through  the  wall  thickness  and  which  has  concrete,  reinforcing  steel 
and  prestressing  steel  layers.  The  contribution  of  each  layer  towards 
the  stiffness  of  the  entire  segment  is  found  individually  and  by  inte¬ 
grating  through  the  volume  of  the  body,  the  stiffness  of  the  segment  as 
a  whole  can  be  found. 

Assuming  a  unit  length,  the  strains  for  each  layer  are  (Bushnell, 

1974) 

en  =  (eimS  -  Z  <h) 

(1  +  Z/Ri) 

e 22  =  (e2mS  ~  z  4>2>  (2.3.11) 

(1  +  Z/R2) 

e  12  = 


0 
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where  Z  is  a  distance  coordinate  measured  from  the  centre  of  the  segment 
thickness,  positive  being  outward  from  the  centres  of  curvature,  and 
Rl  and  R2  are  the  principal  radii  of  curvature  of  the  segment  in  the 
principal  directions  as  shown  in  Fig.  1.1. 

The  element  of  volume  of  one  layer  is 


dV  =  r  (1  +  Z/Ri)  (1  +  Z/R2)  dZ  d0  ds  (2.3.12) 


where  dZ  is  the  thickness  of  the  layer  and  ds  and  rd0  are  the  width  and 
arc  length  of  the  layer,  r  being  the  radius  of  the  parallel  circle  of  the 
reference  surface. 

By  defining  three  new  terms  ki,  k2  and  k3  such  that 


kx  =  (1  +  Z/Rj) 

k2  =  (1  +  Z/R2)  (2.3.13) 

k3  =  klk2 

and  by  assuming  that  e_^  is  associated  only  with  the  mechanical  strains, 

i.e. 


t  c  temp 

e .  =  c .  -  e .  e . 

1111 


(2.3.14) 


t  c 

where  e.  are  the  real  total  strains,  e.  are  the  creep  strains  and 
1  1 

^  temp  a re  the  temperature  strains,  then  differentiation  of  Eq. 


2.3.11  yields 


10 


3e  11 

3e22 

1 

0 

T ''  ms 

de  i 

kl 

3e 

3e  1 1 

3e22 

-z 

0 

3"q 

3<J)  1 

3<j>l 

kl 

3e  1 1 

3e22 

0 

1 

3e2taS 

r\  Ills 
o  £  1 

k  2 

3en 

04 

04 

OJ 

(TO 

0 

-Z 

CM 

-e- 

ex> 

_ 1 

2 

k2 

2 

3  e, 
k 


3q . 3q  . 
i  J 


0 


(2.3.15) 


(2.3.16) 


Substituting  Eqs .  2.3.12,  2.3.13,  2.3.15  and  2.3.16  into  Eq . 
2.3.6  and  carrying  out  the  multiplication 


<Q 


I 


> 


<  ?li 

kl 


Zon  022  ~Zg22  >  kik2  rd0  ds  dz 

ki  k2  k2 


(2.3.17) 


where  <Q  >  are  the  internal  forces  developed  to  resist  the  external 
forces  <Q>.  Assuming  a  unit  length  along  the  reference  surface,  define 
the  size  of  the  segment  as 


rd6  =  1  (2.3.18) 

ds  =  1 


which  results  in 


> 


t 

T 

=  -f  <k20  i  j ,  -Zk20  1 1 ,  k i022 >  ~Zk2<^22>  dZ  (2.3.19) 

-t 

T 


<Nj  Mx  N2  M2> 
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Applying  Eq .  2.3.15  and  2.3.16  to  Eq .  2.3.8,  the  incremental 
equilibrium  expression  becomes 


/ 

1/k ! 

0 

3a  1 1 

3a  u 

3a  ii 

3a  i  i 

dV 

F .  ms— 
Ae  i 

- 

"aNi" 

V 

-Z/k], 

0 

ms 

3e  i 

34>  1 

~  ms 
d£  2 

d<$>2 

A<h 

AMi 

0 

l/k2 

do  22 

do  2  2 

3a22 

3a  2  2 

.  ms 
Ae2 

an2 

0 

-Z/k2_ 

r,  ms 
3e  i 

34>i 

~  ms 
3c  2 

9<J>2 

I 

<N 

-e- 

< 

_ 1 

am2 

0 


(2.3.20) 


Assuming  that 


Oil 

-  f  1 1  (en» 

e  22) 

(2.3.21) 

a22 

=  f 2  2  (£11» 

e  22) 

where  qj  and  £22  are  the  mechanical  strains  as  defined  by  Eq .  2.3.14 
then 


9o 1 1  =  9f  1 1  =  3  f  1 1 


Similarly 


Son  = 

1 

3  f  1 1 

~  ms 

3ei 

kl 

3e  1 1 

9a  1 1  

-Z 

3  f  1 1 

9<h 

kl 

3eil 

9°1 1  

1 

3  f  1  1 

~  ms 

3e  2 

k2 

3e22 

9al  1  

-Z 

3  f  1 1 

3<f>  2 

k2 

9e  22 

9e  11  4. 

3  f 1 1  9a22 

ms 

3ci 

3e  22  9aimS 

3a  22 

1 

3f  22 

in  ms 

3e  1 

kl 

3e  1 1 

3a22 

-Z 

3f  22 

94>1 

kl 

3ei  1 

3a22 

1 

9  f  22 

ms 

3e  2 

k2 

9e  22 

3a  2  2 

-Z 

9f  22 

3<f>2 

k2 

9e  22 

1  ^£n 

ki  3en 


(2.3.22) 
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Substituting  Eqs .  2.3.22  into  2.3.20  results  in  the  following 


C11  C12  C13  C14 

~AN  i“ 

C21  C22  C23  C24 

AM]_ 

C31  C32  C33  C34 

A?2 

AN  2 

C41  C42  C43  C44 

A(J>2 

am2 

[C]  { Aq}  =  {AQ} 


(2.3.23) 


(2. 3. 24) 


where  the  coefficients  of  the  matrix  C  are  defined  as 


t 


C1 1 

7 

=  / 

-t 

7 

9  f  1 1  k  2  dZ 

9en  kl 

C 1  2 

C1  3 

=  / 

9fn  dZ 

9e22 

C 1 4 

C21 

=  f 

— Z  9  f  1 1  k  2  dZ 

9£  1 1  kl 

c22 

c23 

=  / 

-Z  9fn  dZ 

9£  22 

c24 

C  3  1 

=  / 

9  f  2  2  dZ 

9£11 

c32 

C33 

=  / 

9 f 22  k  l  dZ 

3£22  k2 

C34 

C4  1 

=  / 

-Z  9f22  dZ 

9£1 1 

C42 

C43 

=  / 

-Z  9f22  kl  dZ 
9022  k2 

C44 

t 

~Z 

f 

-t 

~1 


9  f  1 1 

k2 

dZ 

9£  1 1 

kl 

9  f  1 1 

dZ 

9£22 

9  f  1 1 

k^ 

dZ 

9  £  1 1 

kl 

9f  1 1 

dZ 

9£22 

3f22 

dZ 

9£1 1 

9f  22 

kl 

dZ 

9£22 

k2 

9f  22 

dZ 

9  £  1 1 

9f  22 

dZ 

9t  22  k2 


(2.3.25) 
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By  solving  Eq.  2.3.24  for  the  unknown  quantities  {Aq}  successive 
corrections  can  be  made  to  {q}  until  {Y}  equals  zero  at  which  point 
the  solution  has  converged. 


2.4  Solution  Technique 

As  stated  in  Section  2.3,  the  Newton-Raphson  iterative  technique 
is  used  to  solve  the  non-linear  coupled  equations  and  follows  the  general 
outline  below. 

1.  Given  an  estimate  of  {q},  calculate  the  forces  {Q^} 
from  Eq .  2.3.19 

2.  Find  the  unbalanced  forces  {AQ}  =  {Q}  -  {Q^} 

3.  If  {AQ}  is  small,  then  the  solution  is  correct.  If  {AQ} 


is  large,  go  to  step  4 

4.  Evaluate  [c]  from  Eqs .  2.3.25 

5.  Find  {Aq}  =  [c]  _1  {AQ} 

6.  New  {q}  =  {q}  +  {Aq} 

7.  Go  to  2 


2.5  Creep  and  Thermal  Loads 

Strains  resulting  from  creep  and  temperature  changes  are  removed 
from  the  total  strains  in  a  layer,  i.e. 

t  c  temp 

e .  =  E  .  £  .  —  £  . 

1111 

t  temp  ,  .  .  . 

where  £.  are  the  total  real  strains,  e.  are  the  strains  resulting 

l  i 

Q 

from  a  temperature  change,  Ei  are  the  strains  in  the  concrete  due  to 
creep  and  e^  are  the  resulting  mechanical  strains.  It  is  the  mechanical 


•: 


• 
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strains  which  are  then  used  to  find  stresses,  tangential  moduli,  and 
Poisson’s  ratios. 


The  variation  in  temperature  changes  through  the  wall  thickness 
can  be  specified  by  one  of  the  following  equations. 


1. 

AT  = 

Tl 

+ 

T2  • 

Y  +  T3  • 

Y2 

2. 

AT  = 

Tl 

+ 

T2  • 

yT  3 

(2.5.1) 

3. 

AT  = 

Tl 

+ 

T2  • 

eY*T3 

Tl, 

T2  and 

T3  are 

user 

defined 

constants 

and  Y  is 

the 

location  coordinate 

The 

origin 

for  Y 

(Y  = 

0) 

is  at 

the  inside 

surface 

of 

the  wall  segment 

(Z  =  -t/2). 

Depending  upon  the  location  and  orientation  of  the  wall  segment 
in  a  structure,  one  or  more  of  the  degrees  of  freedom  of  the  wall  may 
be  restrained  from  temperature  movements.  For  this  reason,  a  program 
variable  must  be  input  to  specify  which  degrees  of  freedom  are  fixed 
against  thermal  movement. 

Creep  strains  are  found  by  using  a  ’long-term'  Young's  modulus. 
During  the  creep  step,  this  reduced  modulus  replaces  the  normal  or 
'short-term'  modulus  and  under  constant  load,  the  strains  increase  to 
maintain  resistance  to  this  load.  The  'long-term'  modulus  is  related 
to  the  'short-term'  modulus  by  the  following. 


g 

where  is  the  creep  factor,  E  is  the  'short-term'  or  normal  Young's 
modulus,  and  EL  is  the  'long-term'  modulus.  The  third  term  applies  for 


■ 


' 

15 


C  0 

a  plain  concrete  section  where  e  would  be  the  creep  strains  and  e 


the  initial  elastic  strains  before  creep. 


The  creep  factor,  F^,  can  be  related  to  the  C.E.B.  creep 
coefficient  (fi  (C.E.B. ,  1969)  by  the  expression 


F 


1 


(2.5.3) 


c 


where 


d>  =  K  K,  K  ,  K  K 
yt  c  b  d  e  t 


(2.5.4) 


and  the  value  is  a  function  of  the  relative  humidity,  is  a  function 

of  the  water  -  cement  ratio,  is  a  function  of  the  age  at  loading, 

K  is  a  function  of  the  fictitious  thickness,  and  K  is  a  function  of 
e  t 


the  duration  of  load.  The  values  of  these  parameters  may  be  obtained 
from  graphs  by  C.E.B.  (1969). 


load.  While  the  program  considers  the  load  on  the  wall  segment  to  be 
kept  constant  during  the  creep  step,  the  load  on  the  concrete  layers 
decreases  as  stress  is  transferred  to  the  steel  layers.  Equation  2.5.2 
therefore  overestimates  the  creep  strain  e^.  This  is  corrected  after 
returning  the  tangential  modulus  back  to  the  short  term  E  by  adjusting 
in  the  concrete  layers  by  successive  iterations  to  compensate  for 
the  reduced  concrete  stress.  The  increase  in  strains  in  the  steel  layers 
which  accompany  the  creep  step  are  real;  the  mechanical  strains  in  the 
steel  layers  equal  the  total  strains. 


3.  CONSTITUTIVE  RELATIONS 


3.1  Introduction 

This  chapter  deals  in  detail  with  the  constitutive  relations  for 
concrete  and  steel.  Two  constitutive  relations  for  concrete  are  pre¬ 
sented.  The  first  is  a  uniaxial  bilinear  elastic  relation  incorporating 
a  tension  cut-off  with  crack  closure  but  no  crack  healing.  Although 
both  directions  on  the  wall  segment  can  be  analysed  using  this  relation, 
they  are  essentially  independent  of  each  other.  The  stress-strain  curve 
of  this  uniaxial  concrete  relation  is  based  on  mechanical  strains.  The 
second  relation  is  a  non-linear  two-dimensional  concrete  constitutive 
relation  which  some  authors  label  hypoelastic  (Schnobrich,  1977).  This 
relationship  is  based  on  the  equivalent  uniaxial  strain  theory  (Darwin 
and  Pecknold,  1974)  and  uses  the  Saenz  equations  (Saenz,  1964)  to  define 
the  stress-strain  law,  in  terms  of  equivalent  uniaxial  strains.  Finally, 
details  of  the  uniaxial  stress-strain  curves  for  both  the  reinforcing 
steel  and  the  prestressing  steel  are  discussed. 

3.2  Uniaxial  Bilinear  Elastic  Concrete  Relation 

The  first  approximation  of  the  concrete  stress-strain  curve  used 

in  this  study  is  a  uniaxial  bilinear  elastic  law  which  requires  the 

specification  of:  E,  Young's  modulus,  f ' ^ ,  the  ultimate  strength  of 

concrete  in  uniaxial  compression,  f ' ^ ,  the  ultimate  strength  of  concrete 

in  uniaxial  tension,  e  ,  ultimate  concrete  strain  in  uniaxial  compres- 

u 

sion,  ec,  which  equals  fWE  and  which  equals  f'^/E.  It  has  been 
shown  (Murray,  1979;  Whitney,  1943)  that  the  post  ultimate  portion  of 
the  concrete  stress-strain  curve  which  can  be  measured  experimentally 
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is  a  function  of  the  stiffness  of  the  testing  machine.  Therefore, 

published  data  concerning  this  property  is  affected  by  the  testing 

procedure.  For  this  model,  it  was  assumed  that  e  =  -  .0038 

u 

(Hognestad  et  al . ,  1955). 

The  following  mathematical  equations  then  describe  the  constitutive 
relation,  which  does  not  account  for  any  Poisson's  ratio  effects. 


a 


a 


a  = 


a 


E  •  e 

0 

f  ' 

c 

0 


E  =  E 


E  = 


E  = 


E  = 


0 

0 

0 


£  <  £  <  £ 
c  t 


£  <  £ 

t 


£  <  £  <  £ 

U  C 


£  <  £ 


U 


(3.2.1) 

(3.2.2) 

(3.2.3) 

(3.2.4) 


This  curve  is  shown  graphically  in  Fig.  3.1.  For  these  and  the  following 
equations,  f?  has  a  value  less  than  zero. 

When  £  exceeds  £  ,  cracking  occurs,  and  for  that  particular  layer, 
the  value  of  z^  is  then  changed  so  that  e  equals  0.  Therefore,  the 
concrete  can  no  longer  develop  tensile  stresses  but  its  behavior  in 
compression  is  assumed  to  be  unaffected.  The  crack  in  this  layer  can 
close  and  compression  stresses  can  form  but  no  subsequent  tensile 
stresses  can  occur  if  the  strains  were  to  become  positive  a  second  time. 

In  compression,  the  stress-strain  curve  is  elastic  up  to  e  . 

Beyond  this  point  however,  the  concrete  in  this  layer  is  assumed  to  have 
crushed  and  can  no  longer  develop  any  stresses  in  tension  or  compression 
regardless  of  the  level  of  strain. 


3.3  Biaxial  Nonlinear  Elastic  Constitutive  Concrete  Relation 

There  are  currently  many  theories  available  for  predicting  the 
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nonlinear  biaxial  behavior  of  plain  concrete  (Chen  and  Chen,  1975; 

Liu  et  al . ,  1972;  Romstad  et  al.,  1974;  Kupfer  and  Gerstle,  1973; 

Bazant  and  Bhat,  1976;  Darwin  and  Pecknold,  1977;  Argyris ,  1973). 

Common  to  most  of  these  models  is  the  assumption  of  linear  elastic 
behavior  of  concrete  in  tension  with  a  tension  cut-off.  Once  cracked, 
the  concrete  is  assumed  to  have  zero  stiffness  in  tension.  In  a  rein¬ 
forced  concrete  specimen  however,  the  cracks  which  form  relieve  the 
tensile  stress  in  the  uncracked  concrete  and  the  concrete  between  the 
crack  continues  to  contribute  some  stiffness  to  the  specimen.  To  account 
for  this  additional  stiffness,  a  few  investigators  have  gone  beyond  the 
simple  tension  cut-off  model  and  have  proposed  a  concrete  stress-strain 
curve  in  tension  complete  with  a  declining  branch.  (Lin  and  Scordelis, 
1975;  Scanlon  and  Murray,  1974;  Murray  et  al,  1978;  Elwi  and  Murray, 

1979)  . 

The  proposed  biaxial  model  incorporates  the  concepts  of  equivalent 
uniaxial  strain  developed  by  Darwin,  Bashur  and  Pecknold  (Darwin  and 
Pecknold,  1974,  1977;  Bashur  and  Darwin,  1978),  the  nonlinear  stress- 
strain  curve  of  Saenz  (1964)  and  the  stress  failure  surface  of  Kupfer 
and  Gerstle  (1973) . 

3.3.1  Equivalent  Uniaxial  Strain 

The  concept  of  equivalent  uniaxial  strains  was  developed  to  remove 
the  Poisson’s  ratio  coupling  of  strains  in  two  (or  three)  dimensions, 
and  therefore  result  in  a  set  of  strains  which  are  independent  of  each 
other,  and  which  can  be  simply  applied  to  a  chosen  stress-strain  law. 

In  this  manner,  it  is  assumed  that  biaxial  stress-strain  curves  can  be 


* 

' 


19 


derived  from  uniaxial  curves.  (Darwin  and  Pecknold,  1974). 

The  incremental  stress-strain  relations  for  an  orthotropic  material 
not  subject  to  shear  take  the  form 


Rail  1 

r Ei 

U2ei" 

“de  1  “ 

|_da  2J  (l-uiu2) 

Lu1E2 

E2  . 

de2 

(3.3.1) 


From  Maxwell’s  reciprocal  theorem  the  stiffness  matrix  should  be  symmetric, 
which  gives  rise  to 


U1E2  -  D2E}  (3.3.2) 

To  simplify  the  expression  further,  the  concept  of  an  ’equivalent 
Poisson’s  ratio'  is  introduced  whereby 

u2  =  u 1U2  (3.3.3) 


Equation  3.3.1  now  takes  the  form 


’dai ' 
da  2 


(1-u2) 


Ei  u  /E 1E2  "|  fde  1' 
u  /E 1 E  2  E  2 


][::;] 


(3.3.4) 


The  incremental  equivalent  uniaxial  strains  are  defined  as 


de  1 


da 


u 


de 


2u 


da; 

E2 


(3.3.5) 


which  upon  substitution  into  Eq .  3.3.4  gives  the  equivalent  uniaxial 
strains  in  terms  of  the  true  strains,  the  elastic  moduli  and  the 
equivalent  Poisson's  ratio. 
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The  total  equivalent  uniaxial  strains  can  now  be  defined  as 

■i«  -  '  i  =  1,  2  (3.3.7) 

E . 

l 

where  the  summation  convention  is  not  used. 

The  equivalent  uniaxial  strains  are  not  real  strains.  They  cannot 
be  observed  and  cannot  be  transformed  under  rotation  of  the  coordinate 
axis  in  the  same  manner  as  stresses  and  true  strains.  Equivalent  uni¬ 
axial  strains  are  accumulated  in  the  principal  directions  which  can 
change  during  loading.  This  theory  does,  however,  recognize  stress- 
induced  orthotropic  behavior.  Fundamentally,  it  allows  the  use  of  curves 
similar  to  that  representing  the  uniaxial  stress-strain  curve  to  describe 
the  biaxial  stress-strain  curves. 

3.3.2  Saenz  Stress-Strain  Curve 

The  uniaxial  stress-strain  relationship  of  Saenz  will  be  extended 
to  cover  biaxial  response,  paralleling  the  procedure  used  by  Darwin  and 
Pecknold  (1977) .  The  same  curve,  with  some  modification,  will  be  used 
to  describe  behavior  in  both  tension  and  compression. 

In  a  uniaxial  test,  real  strains  and  equivalent  uniaxial  strains 
are  equal.  With  this  knowledge,  the  Saenz  equation  will  be  modified 
to  cover  biaxial  response  by  substituting  equivalent  uniaxial  strains 
for  real  strains.  The  Saenz  relation  then  becomes  for  the  direction  i 
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(3.3.8) 


where 


“e 

=  E  /  E 
o  s 

(3.3.9) 
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s 

=  a.  /  e . 
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(3.3.10) 
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(3.3.12) 
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(3.3.13) 

(R  -  l)2  R 

£  £ 

and  where  is  the  current  accumulated  equivalent  uniaxial  strain  in 
direction  i;  a.  is  the  maximum  stress  on  the  stress-strain  curve  in 

1C 

direction  i  for  the  current  principal  stress  ratios;  £  -*-s  t*ie  equiva¬ 

lent  uniaxial  strain  at  a.  ;  o.c  and  £ . _  are  the  coordinates  of  a  point 

ic  if  if 

on  the  declining  branch  of  the  current  stress-equivalent  uniaxial  strain 

curves;  and  E  is  the  initial  modulus  of  elasticity  from  a  uniaxial  test, 
o 

The  type  of  curve  described  by  Eq .  3.3.8  is  shown  in  Fig.  3.2. 

As  required  by  Eq .  3.3.4,  the  incremental  elastic  moduli  must  be 
determined.  From  Eq .  3.3.5  it  can  be  shown  that  for  direction  i 


da . 
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E. 
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(3.3.14) 
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Differentiating  Eq .  3.3.8  with  respect  to  e.  results  in 


da 
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de 


1U 


E. 
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-  2R  /e.  '3 
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(3.3.15) 


1  +  (R+R^-2)  /eiu  \  +  (1-2R)  /e  ^  +  R/e 
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e . 
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Having  determined  the  value  of  E_^ ,  the  change  in  equivalent  uniaxial 
strain,  de .  ,  can  be  calculated  and  accumulated  to  find  the  total 

1U 

equivalent  uniaxial  strain,  e  which  is  required  to  compute  the  stress, 

a .  . 

l 


3.3.3  Poisson's  Ratio 

In  addition  to  the  incremental  elastic  moduli,  Eq .  3.3.4  requires 
the  values  of  Poisson’s  ratio  for  various  locations  on  the  loading  path. 
The  equation  describing  Poisson’s  ratio  was  taken  from  similar  studies 
(Elwi  and  Murray,  1979)  which  used  the  uniaxial  compression  curves  of 
Kupfer,  Hilsdorf  and  Rusch  (1969).  From  this  data,  the  Poisson's  ratio 
was  determined  and  expressed  as  a  function  of  equivalent  uniaxial  strains. 

The  value  for  Poisson’s  ratio,  u.  in  direction  i  is 

l 


U. 

l 


u 

o 


1  +  1.3763 


5.360  /e.  \2  +  8.5  86 

1U 


IUC 


(3.3.16) 


u .  <  0.50 

i 


(3.3.17) 


where  u  is  the  initial  value  of  Poisson’s  ratio.  The  limitation  placed 
o 

on  u.  is  present  to  prevent  the  dilatency  phenomena  observed  in  concrete 
when  reaching  ultimate  strength.  The  equivalent  Poisson’s  ratio,  u. 
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is  now  found  from  Eq .  3.3.3. 

3.3.4  Ultimate  Curves 

To  define  the  stress-strain  curve  described  by  Eq.  3.3.8,  the 

additional  parameters  accompanying  this  equation  must  be  determined. 

These  parameters  are  a.  ,  a.-,  e.  and  e.  _.  a.  and  e.  will  be 

ic  if  iuc  luf  1C  1UC 

expressed  as  failure  curves  which  are  direct  or  indirect  functions  of 

the  current  principal  stress  ratio,  and  o . c  and  e._  will  be  functions 

if  if 

of  a.  and  e. 

1C  IUC 


3. 3.4.1  Ultimate  Stress  Curve 

The  stress  failure  curve  of  Kupfer  and  Gerstle  (1973)  was  used  to 
determine  the  values  of  a  .  This  curve  is  composed  of  three  parts 
which  are  functions  of  the  stress  ratio  and  is  illustrated  in  Fig.  3.3. 

In  the  compression-compression  quadrant  of  the  biaxial  stress 
space,  the  strength  envelop  is  described  by 

(°lc  +  a2cV  -  a2c  -  3.65  °lc  =  0  (3.3.18) 

If’  f '  1  f f  f ' 

\  C  cj  c  c 

for  which  is  greater  than  or  equal  to  a 2, and  o^,  o 2  and  f*  are  values 
less  than  zero. 

By  defining  the  ratio  as 

°J_  (3.3.19) 

Eq.  3.3.18  can  be  solved  to  express  Oj  and  a2c  in  terms  of  as 
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1  +  3.65  a  £ , 

(i  +  «)2 


(3.3.20) 


a 


1  c 


(3.3.21) 


In  this  case,  a2c  is  the  major  direction  (the  direction  of  higher  absolute 
stress) . 

In  the  tension-compression  quadrants,  the  following  straight  line 
is  used. 
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lc 


f 
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1  -  0.8  o0 

zc 

f  ’ 
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which  can  be  restated  as 


«f'  +0 . 8  f ’ 
c  t 


a 


1  c 


ac 


a,  >  0 

1  c 


(3.3.22) 


(3.3.23) 


(3.3.24) 


There  is  a  small  region  of  overlap  where  Eqs .  3.3.20  and  3.3.21  are  used 
in  the  tension-compression  quadrants. 

For  the  tension-tension  quadrant,  a  constant  tensile  strength  is 

used . 


a,  =  f’  (3.3.25) 
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(3.3.26) 
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3. 3. 4. 2  Equivalent  Uniaxial  Strain  at  Ultimate  Stress 

The  equations  which  were  derived  for  this  section  are  based  on 
the  data  from  the  biaxial  tests  of  Kupfer  (1973) .  A  program  developed 
by  A.  A.  Elwi,  using  the  stress-strain  relations  for  an  orthotropic 
material,  the  Saenz  (1964)  stress-strain  equations  and  the  Poisson’s 
ratio  equation  of  A.  A.  Elwi,  has  been  used  to  mathematically  model 
the  constant  stress  ratio  experiments  of  Kupfer.  By  varying  the  input 
value  of  e .  into  this  program,  the  optimum  e .  which  provided  a 

1UC  1UC 

reasonable  correlation  between  the  predicted  and  observed  stress-strain 

curves  was  found  for  each  stress  ratio.  Figures  3.4  through  3.8  compare 

the  curves  of  Kupfer,  for  a  few  stress  ratios,  with  those  obtained  using 

the  above  technique  and  the  optimum  values  of  e. 

iuc  . 

The  values  of  e.  determined  to  optimize  the  correlation  with  the 

iuc 

specific  set  of  stress-strain  curves  provided  by  Kupfer  (1973) ,  were 
then  generalized  to  define  an  equivalent  uniaxial  strain  curve  in  biaxial 
space.  That  is,  from  the  expressions  of  e  which  were  found,  expres¬ 
sions  were  developed  to  relate  the  current  stress  ratio  to  the  associated 

value  of  e.  ,  where  e.  is  the  equivalent  uniaxial  strain  when  the 
iuc  iuc 

stress,  o^  equals  the  maximum,  cj^  ,  for  that  particular  stress  ratio. 
Figure  3.9  contains  a  diagram  of  the  e  envelop. 

In  the  compression-compression  region,  a  new  variable  6  was  defined 

where 

6  =  £luc  (3.3.27) 


A  new  material  constant  e_^o,  is  now  required  to  define  the  curve 
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in  Fig.  3.9,  where  e  is  the  equivalent  uniaxial  strain  in  the  stressed 
direction  at  ultimate  stress  for  a  uniaxial  compression  test.  e.  can 

1UC 


now  be  calculated  using 


£  2uc 

.  elo  (1  +  7.756) 

(3.3.28) 

(1  +  46) 

£  luc 

=  3 

2uc 

(3.3.29) 

6 

=  cc3  -  1.8  cc2  +  1.8  cc 

(3.3.30) 

When  02  equals  f ' 

and  oi  equals  zero,  then  equals  e.  . 

A  second  material  constant,  e.  is  required  where  e.  equals  the 

lot  lot 

equivalent  uniaxial  strain,  in  the  stressed  direction,  at  ultimate  stress, 
for  a  uniaxial  tension  test.  In  the  tension-compression  quadrants,  the 
following ‘relationship  is  suggested. 


When  Oi  equals  f'  and  a?  equals  zero,  e,  =  e. 

1  M  t  Z  M  1UC  lot 

For  the  tension-tension  quadrant,  e.  is  found  by 

1UC 
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1  uc 
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2uc 


£iot 
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(1  +  0.5  0  -  0.25  et2) 

(0.75  4-  0.5  Bt3) 
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-  1  /  <x 


(3.3.31) 

(3.3.32) 


(3.3.33) 

(3.3.34) 

(3.3.35) 
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Table  3.1  gives  a  comparison  of  the  observed  values  of  a.  and 

1C 

e .  against  calculated  values  for  a  medium  strength  concrete.  In  this 

1UC 

table,  the  stress  and  strain  values  are  nondimensionalized  with  the 

values  of  f1  ,  f’  ,  e.  and  e.  .  ,  which  for  this  concrete  were  -4650  psi, 
c  t  10  lot  r 

423  psi,  -.00219  in. /in.  and  .00009  in. /in.  respectively.  The  observed 

values  are  those  taken  from  Kupfer  (1973) .  The  calculated  values  of 

a.  /  f*  and  a  /  f*  are  the  values  which  would  be  found  using  the 
1C  c  1C  t 

equations  in  Section  3. 3. 4.1  and  the  given  stress  ratios.  Two  sets  of 

the  values  e.  /  e.  and  e.  /  e.  are  then  listed.  The  first  set, 
iuc  10  iuc  lot 

calculated  from  the  predicted  a  represent  the  values  which  would  be 
found  starting  only  with  the  knowledge  of  the  stress  ratio,  and  include 
any  errors  which  may  be  included  when  calculating  .  The  second  set 
are  calculated  using  the  observed  values  of  a.  .  However,  there  is  no 

1C 

difference  in  the  strain  values  obtained  using  the  predicted  or  observed 
0L  in  the  compression-compression  or  tension-tension  quadrants  as  the 
strain  values  in  these  quadrants  are  calculated  directly  from  the  stress 
ratio  and  independently  of  the  ultimate  stresses.  Only  in  the  compres¬ 
sion-tension  region  are  the  ultimate  stresses  used  to  calculate  e. 

&  iuc 

The  values  of  a.f  and  are  used  to  define  the  shape  of  the 

if  if 

declining  branch  of  the  stress-strain  curve.  At  the  present  time,  there 
are  few  existing  studies  on  the  declining  branch  of  concrete  stress- 
strain  curves  (Evans  and  Morathe,  1968;  Hughes  and  Chapman,  1966)  and 
the  following  assumptions,  made  by  Elwi  and  Murray  (1979)  will  be  used. 


e .  £  =  4  e . 

if  iuc 

a.  _  =  a.  /  4 

if  1C 


(3.3.36) 


(3.3.37) 
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These  are  the  last  values  required  to  calculate  and  E_^  from  Eqs . 

3.3.8  and  3.3.15. 

3.3.5  Tensile  Declining  Branch 

Chapter  5  contains  a  detailed  comparison  of  the  experimental 
observations  of  test  wall  segments  and  of  the  analytical  results  obtained 
from  COWSAC  on  models  of  these  test  segments.  From  these  comparisons, 
it  became  clear  that  the  declining  branch  of  the  Saenz  stress-strain 
curve  was  too  steep  for  tensile  stresses.  As  a  result,  the  declining 
branch  previously  described  by  the  Saenz  equation  was  replaced  by  a 
straight-line  of  negative  slope  with  a  minimum  residual  strength.  This 
curve  is  shown  graphically  in  Fig.  3.10.  Up  to  e  ,  in  both  the  major 
and  minor  directions,  the  a  -  e.  curves  follow  the  Saenz  relation. 

1U 

Beyond  these  points,  the  stress  decrease  Ack,  is  linearly  related  to  the 
increase  in  equivalent  uniaxial  strain,  Ae .  . 

1U 

For  the  major  direction,  (direction  1,  0]_  >  a 2)  the  declining  slope 

Ei  is 


Ei  =  6.5a 


1, 


7.5  x  30  e . 

1UC 


In  the  minor  direction,  the  slope  E2  equals 


(3.3.38) 


E9  =  6 .5  a9 

^  _ £C 

7.5  x  30  e. 

1UC 


(3.3.39) 


In  both  directions  the  minimum  residual  stress  is  a,  /  7.5. 

1  c 
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As  the  Saenz  equation  was  developed  for  concrete  subject  to  compression 
forces,  and  no  test  wall  segments  were  tested  in  axial  compression,  there 
was  no  necessity  to  change  the  postultimate  portion  of  the  Saenz  curve 
in  compression.  For  tension-compression  stress  combinations,  the  tensile 
stresses  are  assumed  to  follow  the  a  -  e.  curve  of  the  major  direction 

1U 

shown  in  Fig.  3.10  while  the  compressive  stresses  follow  the  Saenz 
equation. 


3.3.6  Equivalent  Uniaxial  Strain  Constants 

To  accurately  describe  the  concrete  constitutive  relations,  two 

additional  constants  over  and  above  the  normally  obtainable  values  of 

ff  ,  f’t,  Eq  and  l>o  are  required  to  help  shape  the  concrete  stress-strain 

curves  in  compression  and  tension.  These  are  e.  ,  the  strain  at  maximum 
r  10 

stress  in  uniaxial  compression,  and  e.  ,  the  strain  at  maximum  stress  in 

lot 

uniaxial  tension.  For  uniaxial  tension  or  compression,  the  equivalent 
uniaxial  strains  equal  the  real  strains. 

If  these  strains  cannot  be  provided  to  the  program,  the  following 
default  equations  are  invoked  (Saenz,  1964) 


e  . 
10 


-1  x  f’ 


(31.5 


100,000 


(3.3.40) 


Tot  (3.3.41) 

E 

o 

For  Eq.  3.3.40,  f’  must  be  in  units  of  psi.  The  value  of  e.  from 

1  C  IOC 

Eq.  3.3.41  assumes  that  the  a  -  eiu  curve  in  tension  is  linear  up  to  the 
ultimate  stress  whereas  the  program  follows  the  curved  Saenz  equation, 
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whose  slope  starts  at  E  from  a  strain  of  zero  and  then  declines  to 

o 

zero  at  a  strain  equaling  .  Therefore,  it  is  suggested  that  the 

default  value  of  e.  be  overridden  when  using  COWSAC.  A  better  value 

lot 

of  e  ^  would  be  larger  than  that  suggested  by  Eq .  3.3.41,  possibly  by 
a  factor  of  1.2  or  1.3. 

3.3.7  Instability  Arising  From  the  Declining  Branch  of  the  Stress-Strain 

Curve 

The  use  of  the  concept  of  equivalent  uniaxial  strains,  and  of  the 
stress  ratio,  aj  /  a 2  as  the  basis  for  defining  the  failure  curves,  in 
conjunction  with  an  iterative  numerical  procedure,  leads  to  instabilities 
on  the  declining  branches  of  the  a  -  e.  curves.  This  instability 

1U 

manifests  itself  in  the  degeneration  of  the  minor  stress-strain  curve 
(smaller  absolute  magnitude)  to  a  point  at  the  origin  of  the  stress- 
equivalent  uniaxial  strain  space. 

To  describe  this  problem,  some  of  the  solution  procedures  will  be 
described  for  test  wall  segment  No.  1,  detailed  in  Table  5.1  and  Fig.  5.4. 
This  is  a  typical  prestressed  and  reinforced  wall  segment  tested  at  a 
load  ratio  of  2:1.  For  this  specimen,  the  direction  with  the  greater 
load  will  be  called  direction  1.  This  direction  also  corresponds  with 
the  major  direction  (direction  of  greater  absolute  stress)  of  the 
concrete  layers. 

As  indicated  by  the  test  results  in  Fig.  5.4,  each  of  the  two 
load-extension  curves  can  be  divided  into  two  distinct  parts:  the  first, 
of  large  constant  stiffness  up  to  cracking,  in  which  both  the  stress  in 
the  steel  and  the  concrete  increase,  and  the  second  of  smaller  stiffness 
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(slope)  after  cracking,  in  which  the  stress  in  the  concrete  decreases. 

On  the  first  portion  of  these  curves,  when  the  total  specimen 
displays  a  large  stiffness,  the  concrete  is  on  the  ascending  branch  and 
the  minor  direction  is  stable.  If  a  greater  stress  is  required  in 
direction  2,  the  strains  in  that  direction  are  increased.  When  this 
occurs,  the  stress  point  of  the  concrete  in  an  individual  layer,  which 
is  originally  at  location  B  on  the  solid  curve  of  Fig.  3.11a  moves  up 
and  along  this  solid  curve  to  increase  the  stress  in  direction  2.  Since 
this  stress  has  increased  relative  to  the  major  direction,  the  stress 
ratio  changes  and  the  final  location  of  the  stress-strain  point  in 
direction  2  is  at  B'  on  the  dashed  line.  Since  the  total  strains  have 
increased,  the  stress  in  the  steel  layers  also  increases.  For  the  entire 
segment,  the  response  to  an  increase  in  load  in  direction  2  is,  due  to 
the  positive  stiffness  of  the  segment,  an  increase  in  strain  in  direction 
2.  This  results  in  greater  stress  in  the  steel  layers,  and  greater  stress 
in  the  concrete  as  the  concrete  stress  moves  from  B  to  B 1  in  Fig.  3.11a. 

On  the  second  branch  of  the  curves  in  Fig.  5.4,  when  the  total 
specimen  displays  a  smaller  positive  stiffness,  the  concrete  is  on  the 
declining  branch  of  the  stress-strain  curves  and  the  minor  direction 
is  unstable.  If  a  greater  stress  is  required  in  direction  2,  the  strains 
in  direction  2  are  increased,  as  the  overall  segment  has  a  positive  stiff¬ 
ness.  When  this  occurs,  the  stress-strain  point  of  the  concrete  in  an 
individual  layer  which  is  originally  at  location  B  of  Fig.  3.11b  moves 
down  along  the  stress-strain  curve.  Since  the  stress-ratio  has  changed, 
the  final  location  of  the  stress-strain  point  in  the  minor  direction  is 


i  ■ 

' 


32 


at  B'  on  the  dashed  curve  of  Fig.  3.11b.  The  overall  segment  response 
indicates  an  increase  in  steel  stresses,  but  a  large  decrease  in  concrete 
stresses.  To  compensate  for  the  decrease  in  concrete  stresses,  the  strains 
in  direction  2  are  increased  further.  The  cycle  is  repeated  and  the  minor 
direction  degenerates  to  zero. 

To  prevent  this  degeneration,  the  stress  loss  in  absolute  terms 
in  the  minor  direction  on  the  descending  branch  has  been  limited  to 

I  A  ^  I  I  t-i^^  A  ^  I  /  O  O  /  O  \ 

I Aa  |  <  |  E  x  Ae  |  (3.3.42) 

where  Aam  is  the  change  in  stress  in  the  minor  direction,  E™  is  the 
tangent  modulus  on  the  declining  branch  in  the  minor  direction,  and 
Aem  is  the  change  in  real  mechanical  strains  in  the  minor  direction. 

If  the  stress  loss  calculated  by  the  normal  procedure  is  greater 
than  that  calculated  by  Eq .  3.3.42,  then  the  results  are  overridden 
by  Eq .  3.3.42.  This  action  prevents  the  stress  ratio  on  the  declining 
branch  from  decreasing,  but  will  not  inhibit  an  increase  in  the  stress 
ratio,  which  is  updated  twice  during  every  increment  in  real  strains. 

(See  Chapter  4) . 

3.3.8  Unloading 

In  a  generally  accepted  definition  of  unloading  as  used  in  plasticity 
(Fung,  1965)  a  loading  surface  exists  which  is  a  function  of  the  stress 
state,  plastic  strains  and  a  strain  hardening  parameter,  on  which  the 
current  biaxial  stress  state  is  located.  It  can  be  assumed  that  the 
shape  of  this  curve  in  two  dimensions  before  plastic  strains  have  occurred 
is  similar  to  that  describing  the  stress  failure  curve,  and  that  if  the 
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distance  from  the  point  of  zero  stress  to  the  failure  curve  in  any  one 
direction  is  L^,  and  if  the  distance  from  the  point  of  zero  stress  to  the 
loading  curve  in  this  same  direction  is  L^,  then  the  ratio  is  a 

constant  (less  than  or  equal  to  one)  for  any  particular  loading  curve. 
Knowing  the  shape  of  the  loading  curve  then,  in  the  elastic  range  if  the 
next  stress  point  is  on  the  outside  of  the  loading  curve,  the  material 
is  being  loaded.  If  the  next  stress  point  is  on  the  inside  of  the  load¬ 
ing  curve,  then  the  material  is  being  unloaded. 

To  reduce  the  complexity  of  applying  this  definition  in  programming, 
the  failure  curve  used  was  the  curve  of  equivalent  uniaxial  strains  at 
ultimate  stress,  Fig.  3.9,  and  the  loading  curve  contained  the  current 
equivalent  uniaxial  stresses.  Unloading  again  occurs  when  the  next 
point  is  within  the  current  loading  curve.  The  unloading  condition  is 
then  approximated  by  the  following. 


£  . 

e . 

1U 

— 

1U 

<  0 

£  . 

1UC 

present  load  state 

£  . 

1UC 

previous  load  state 

(3.3.43) 

If  Eq .  3.3.43  is  satisfied  in  both  directions,  then  an  unloading  condition 
is  assumed  to  have  occurred.  If  the  material  is  on  the  ascending  portion 
of  the  stress-strain  curve,  then  unloading  follows  the  nonlinear  Saenz 
relation.  However,  if  unloading  occurs  when  on  the  descending  branch 
of  the  stress-strain  curve,  then  the  stresses  and  stiffness  drop  to 
zero.  When  this  happens  on  a  compression  curve,  that  layer  of  concrete 
is  assumed  to  have  crushed  and  will  take  no  further  compressive  or  tensile 
loads.  When  this  occurs  on  a  tensile  curve,  that  layer  of  concrete  is 
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assumed  to  have  cracked  and  will  take  no  further  tensile  stresses.  This 
cracked  layer  can  take  full  compressive  load  after  the  crack  closes.  The 
example  problem  discussed  in  Chapter  4  discusses  unloading  behavior  in 
further  detail. 

3.3.9  Crack  Closure  After  Unloading 

In  a  reinforced  concrete  wall  segment  subjected  to  tensile  loads, 
the  average  response  of  the  concrete  appears  to  follow  the  stress-strain 
curve  exhibited  in  Fig.  3.10.  However  cracks  open  up  in  the  specimen, 
across  which  no  stresses  can  exist,  even  though  the  'average'  stress- 
strain  curve  predicts  some  load  carrying  capacity.  It  can  be  assumed 
that  the  classical  brittle  linear  elastic  curve  with  a  tension  cut-off 
shown  in  Fig.  3.14  more  closely  represents  the  concrete  behavior  on  a 
microscopic  level  in  a  normal  test  situation.  With  this  model,  no  in¬ 
elastic  strains  will  accumulate  in  the  concrete  even  though  the  'average' 
strains  are  greater  than  the  strain  at  maximum  stress.  The  additional 
strains  predicted  by  Fig.  3.10  result  from  the  opening  and  widening  of 
cracks.  The  concrete  in  between  the  cracks  behaves  linearly  and 
elastically.  Therefore  we  can  predict  the  point  at  which  the  concrete 
accepts  compressive  stresses  (crack  closure)  with  elastic  theory,  as 
follows . 

If  the  concrete  is  cracked  in  both  directions,  then  the  concrete 
will  close  in  one  direction  when 

e.  <  0  (3.3.44) 

l 

where  e.  are  the  real  mechanical  strains  in  direction  i.  This  then 
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reduces  the  problem  to  one  of  cracking  in  one  direction. 

If  the  concrete  is  cracked  in  only  one  direction,  the  point  at  which 
the  cracks  close  can  be  found  by  determining  the  conditions  under  which 
the  concrete  will  develop  compression  in  this  direction. 

If  direction  1  is  the  cracked  direction,  then  at  crack  closure, 

Eq .  3.3.4  can  be  evaluated  by  assuming  that  E2>  the  tangent  modulus  in 
the  stressed  direction,  is  approximately  equal  to  the  secant  modulus, 
and  hence 


Ei  =  E 

o 

E  2  =  E2 

Then  Eq .  3.3.4  becomes 


Oi  =  0  (3.3.45) 

02  =  o2  (3.3.46) 


<*l 


(1-u2) 


^  El  £]_  +  U  /  E1E2  £2 


(3.3.47) 


Setting  0]_  to  zero  in  Eq .  3.3.47,  the  strain,  e  1 ,  at  which  cracks  in 
direction  1  close  can  be  determined  as 


=  -  u  /  E2  /  Eq  £2  (3.3.48) 

If  the  second  direction  is  in  compression  and  has  progressed  sig¬ 
nificantly  along  the  stress-strain  curve,  some  inelastic  strains  will 
have  occurred,  violating  the  assumption  of  elastic  behavior  and  reducing 


the  accuracy  of  this  relation. 
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3.4  Reinforcing  Steel  Constitutive  Relation 

The  reinforcing  steel  is  modelled  as  a  uniaxial  layer,  contributing 
its  stiffness  to  the  wall  segment  only  along  its  length,  and  replacing 
the  concrete  in  that  direction.  In  the  transverse  direction,  the  presence 
of  the  steel  has  no  effect  and  the  concrete  stiffness  and  stresses  are 
used,  making  no  modification  for  the  steel.  Due  to  its  one-dimensional 
contribution,  a  simple  elastic-plastic  stress-strain  curve  can  be  used 
for  the  reinforcing  steel  as  shown  in  Fig.  3.12. 

The  steel  has  an  initial  constant  stiffness  up  to  yield,  beyond 
which  the  tangent  stiffness  becomes  zero  and  the  stress  is  maintained 
at  the  yield  stress.  This  is  described  by  the  equations 


a 


£  < 


(3.4.1) 


a  =  E  _  x  e 
st 


<  £ 


(3.4.2) 


a 


a 

y 


£  <  £ 
y 


(3.4.3) 


where  a  is  the  yield  stress  of  the  steel,  E  is  the  Young’s  modulus  for 
y  s  t 

the  steel  and  e  equals  a  /  E  This  curve  is  followed  for  both  loading 

y  y  st 

and  unloading. 


3.5  Prestressing  Steel  Constitutive  Relation 

The  prestressing  strands  are  also  one-dimensional  elements  and  are 
treated  similarly  to  the  reinforcing  steel  elements.  The  stress-strain 
curve  for  the  strands  was  derived  from  tests  taken  in  conjunction  with 
the  wall  segment  tests  discussed  in  Chapter  5  and  is  composed  of  6  straight 
lines  as  shown  in  Fig.  3.13. 

The  equations  describing  the  curve,  again  applicable  to  both  loading 
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and  unloading,  are 

a  =  E  x  e 
P 

a  =  0.777  a  +  (.087  a  )  (e  -  e 

pu  pu 

.0084-e 


y 

a  =  0.864  a  +  20  a  (e  -  .0084) 

pu  pu 

a  =  0.896  a  +  6.5  a  (e  -  .010) 

pu  pu 

a  =  0.909  a  +  4.75  a  (e  -  .012) 

pu  pu 

a  =  0.947  a  +  .053  a  (e  -  .020) 

pu  pu 


.23 


a  =  0 


£  <  e 

sy 


£  <  e  <  .0084 

sy 


.0084  <  £  <  .010 


.010  <  £  <  .012 


.012  <  £  <  .020 


.020  <  £  <  .25 


.25  <  £ 


(3.5.1) 

(3.5.2) 

(3.5.3) 

(3.5.4) 

(3.5.5) 

(3.5.6) 

(3.5.7) 


where  is  the  initial  Young’s  modulus  of  the  prestressing  strands, 

a  is  the  ultimate  strength  of  the  prestressing  strands  and  £  equals 
pu  sy 


0.777  a  /  E  . 
pu  p 


. 
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Observed  Values  Calculated  Values 


From  Pre-  From  Ob- 
dicted  a.  served  a 


1C 

1C 

<*1 

o. 

1C 

a. 

1C 

£  . 

1UC 

E  . 

1UC 

a . 

1C 

a . 

1C 

£  . 

1UC 

Ciuc 

£  . 

1UC 

£  . 
1UC 

02 

f  ’ 

C 

fT 

t 

E  . 

10 

eiot 

f’ 

C 

f ' 

t 

e . 

lO 

eiot 

£  . 
lO 

Eiot 

-1 

1.15 

1.7  35 

1.16 

1.750 

1.750 

-1 

1.15 

1.735 

1.16 

1.75 

1.750 

-0.52 

0.645 

0.890 

0.652 

0.978 

0.978 

-1 

1.240 

1.667 

1.25 

1.658 

1.658 

CM 

CN 

• 

O 

1 

0.262 

1.594 

0.269 

0.488 

0.488 

-1 

1.19 

1.575 

1.21 

1.526 

1.5  26 

0 

0. 

0 

0 

0 

0 

-1 

1.00 

1.00 

1.00 

1.00 

1.00 

0.052 

0.484 

1.11 

0.417 

0.884 

0.929 

-1 

0.850 

0.653 

0.729 

0.493 

0.675 

0.070 

0.594 

1.11 

0.490 

0.932 

0.973 

-1 

0.770 

0.557 

0.637 

0.39  3 

0.548 

0.103 

0.704 

0.667 

0.5  86 

0.971 

0.992 

-1 

0.620 

0.388 

0.518 

0.299 

0.377 

0.202 

0.824 

0.773 

0.735 

0.995 

0.999 

-1 

0.370 

0.210 

0.331 

0.199 

0.217 

1 

1.00 

1.00 

1.00 

1.00 

1.00 

0 

0 

0 

0 

0 

0 

1 

1.077 

1.17 

1.00 

1.10 

1.10 

0.23 

0.253 

0.556 

0.230 

0.466 

0.466 

1 

1.033 

1.11 

1.00 

1.20 

1.20 

0.54 

0.561 

0.667 

0.540 

0.6  89 

0.6  89 

1 

1.00 

1.22 

1.00 

1.25 

1.25 

1 

1.00 

1.22 

1.00 

1.25 

1.25 

TABLE  3.1 


Comparison  of  Observed  and  Calculated  Values  of  a.  and  e. 
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Fig.  3.1 


Uniaxial  Bilinear  Elastic  Stress-Strain  Law 
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Fig.  3.2 


Typical  Compressive  Stress-Equivalent  Uniaxial  Strain  Curve 
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Fig.  3.3 


Ultimate  Stress  Curve 
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Comparison  of  Saenz  Equation  and  KHR  Data  for  0^/02 
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Fig.  3.7  Comparison  of  Saenz  Equation  and  KUR  Data  for  Oj 
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Fig.  3.8  Comparison  of  Saenz  Equation  and  KHR  Data  for  oj 
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Fig.  3.9 


Equivalent  Uniaxial  Strains  at  Ultimate  Stress 
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Fig.  3  .10  Stress-Equivalent  Uniaxial  Strain  Curve  in  Tension 
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Fig.  3.11a  Concrete  Stresses  on  Ascending  Branch  of  a-e^  Curve 


Fig.  3.11b  Concrete  Stresses  on  Descending  Branch  of  a-eiu  Curve 
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Fig.  3.12  Stress-Strain  Curve  for  Reinforcing  Steel 


Stress 


51 


Fig.  3.13  Stress-Strain  Curve  for  Prestressing  Strands 
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Fig.  3.14  Linear  Elastic  Stress-Strain  Curve  for  Concrete  in 

Tension,  with  Tension  Cut-off 


4 .  PROGRAM  DEVELOPMENT 


4.1  Description  of  Approach 

The  computer  coding  which  was  developed  to  solve  the  range  of  pro¬ 
blems  outlined  in  Chapter  1,  was  designed  to  be  user  oriented.  To 
satisfy  this  requirement,  the  successful  operation  of  the  program  must 
be  accomplished  with  a  minimum  of  input,  and  the  data  requests  must  be 
sufficiently  flexible  to  handle  a  wide  variety  of  problems.  To  help 
provide  the  needed  versatility,  no  restrictions  are  placed  on  the  specimen 
dimensions  or  steel  layer  locations.  A  variety  of  loading  schemes  are 
available,  which  can  be  mixed.  Any  consistent  system  of  units  may  be 
used  to  specify  the  problem;  the  program  output  will  follow  in  consistent 
units . 

Before  the  live  load  is  applied,  there  is  opportunity  to  place  four 
different  types  of  load  on  the  wall  segments  in  the  following  order, 
which  is  based  on  a  possible  construction  sequence.  Before  the  pre¬ 
stressing  is  applied,  a  dead  load  can  be  placed  on  the  segment.  This 
load  will  stress  the  concrete  and  reinforcing  steel,  but  not  the  pre¬ 
stressing  steel.  For  this  load  step,  the  area  occupied  by  the  prestres¬ 
sing  steel  will  be  assumed  to  have  no  stiffness  and  no  strength,  so 
effectively  the  prestressing  strands  are  considered  to  be  absent,  with 
empty  holes  at  their  designated  locations.  For  this  and  all  load  steps, 
the  concrete  strength  is  expected  to  have  the  set  of  material  properties 
which  are  specified  by  the  user.  No  reduction  in  strength  due  to  con¬ 
crete  immaturity  is  here  allowed  or  accounted  for. 

After  the  dead  load,  the  prestressing  can  be  applied.  The  pre- 
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stressing  force  (rather  than  stress)  can  be  specified  as  a  value  to  be 
achieved  either  before  or  after  the  elastic  response  of  the  specimen  to 
the  prestressing  forces  has  taken  place.  In  addition,  the  option  exists 
to  sequence  the  prestressing.  The  tendons  can  be  stressed  either 
individually,  in  groups,  or  all  at  once,  according  to  the  order  specified 
by  the  user.  As  each  of  the  prestress  steps  is  applied,  a  load  must  be 
specified  which  is  the  load  on  the  segment  at  the  conclusion  of  each 
prestress  increment.  This  load  is  in  addition  to  the  internal  load 
resulting  from  the  prestressing.  For  example,  the  test  wall  segments 
which  are  discussed  in  Chapter  5  are  prestressed  but  are  not  subjected 
to  any  external  forces.  Therefore,  the  load  on  the  segment  at  the  con¬ 
clusion  of  prestressing  of  these  test  wall  segments  would  be  zero  forces 
and  zero  moments.  However,  the  wall  segments  discussed  in  Chapter  6, 
which  form  part  of  the  containment  structure,  have  gravity  loads  upon 
them.  The  segment  must  be  in  equilibrium  with  the  prestressing  and 
these  specified  gravity  loads. 

Once  the  dead  load  and  prestressing  effects  have  been  determined, 
the  specimen  can  be  allowed  to  creep  if  desired.  This  is  accomplished 
by  changing  the  elastic  modulus  of  the  concrete  from  the  original  short¬ 
term  modulus  to  the  long-term  values,  by  means  of  a  creep  factor  which 
is  discussed  in  Chapter  3.  The  loads  on  the  specimen  at  the  time  of 
prestressing  are  maintained  and  the  creep  strains  which  accumulate  in 
the  concrete  are  recorded.  After  the  creep  step,  the  elastic  modulus 
of  the  concrete  returns  to  its  short-term  value  and  the  creep  strains 
of  each  concrete  layer  are  subtracted  from  the  total  strains  to  result 


in  the  net  mechanical  strains.  The  additional  strains  in  the  steel 
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layer  which  result  from  the  creep  process  remain  as  the  total  strains 
and  are  used  to  calculate  the  steel  stresses. 

The  first  type  of  live  load  which  can  now  be  applied  is  a  thermal 
load.  Three  different  types  of  thermal  gradients,  each  being  defined 
by  three  specified  constants  as  indicated  in  Eqs .  2.5.1  are  available  to 
describe  most  expected  temperature  distributions.  The  temperature  change 
(from  zero)  calculated  at  each  layer  is  then  multiplied  by  the  appropriate 
thermal  expansion  coefficient  to  find  the  thermal  strain  for  that  layer. 
The  thermal  effects  can  occur  under  boundary  conditions  of  constant  load 
or  constant  strain,  as  specified  by  the  user. 

The  final  loads  to  be  applied  are  the  live  loads.  These  loads  can 
be  specified  either  as  the  total  loads  on  the  segment,  or  as  additional 
loads  to  be  added  to  the  previous  total  loads.  The  type  of  load, 
either  total  or  incremental,  can  be  changed  for  each  live  load.  Each 
live  load  can  be  subdivided  into  load  increments  so  that  a  solution 

will  be  found  at  intervals  between  loads.  The  number  of  these  sub¬ 
divisions  can  be  changed  for  each  live  load.  If  so  desired,  a  loading 

procedure  can  be  invoked  where  the  load  consists  of  axial  forces  and 

curvatures,  the  two  being  linearly  dependent,  and  the  curvature  being 
calculated  from  initial  force  and  moment  data.  This  is  the  type  of 
loading  employed  by  Epstein  and  Murray  (1976). 

The  above  loading  sequence  of  dead  loads,  prestressing  loads, 
creep  loads,  thermal  loads  and  live  loads  cannot  be  altered.  However, 


not  all  of  the  loads  must  be  included. 


. 
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4.2  Program  Flowcharts 

The  flowcharts  in  Figs.  4.1  to  4.5  serve  to  provide  a  general 
description  of  logical  flow  within  the  program,  and  to  disclose  the 
main  iterative  loops  and  subincrement  generators.  A  detailed  User’s 
Manual,  specifying  input  requirements,  is  contained  in  Appendix  A,  and 
a  program  listing  is  contained  in  Appendix  B. 

The  main  routine  calls  for  the  solution  of  the  loads  described 
above  in  the  stated  sequence.  This  is  illustrated  in  Fig.  4.1.  The 
solution  of  each  type  of  load  is  similar  and  can  be  traced  through  the 
flowcharts  of  load  control  and  load  application,  Fig.  4.2  and  4.3. 

The  flowchart  of  load  control,  shown  in  Fig.  4.1,  describes  the 
prime  generator  of  load  subincrements.  If  the  bilinear  elastic  concrete 
constitutive  relation  is  used  for  the  solution,  the  load  is  not  sub¬ 
divided.  If  the  nonlinear  concrete  constitutive  equations  are  specified 
however,  then  the  load  is  subdivided  to  restrict  the  change  in  overall 
segment  strains  to  a  set  value.  Once  the  load  is  subdivided,  all  the 
values  describing  the  current  status  of  the  material  layers  are  noted. 
These  values  are  known  to  be  correct  and  shall  be  used  until  the  solution 
has  converged  at  the  next  load  subincrement.  The  method  of  solution  is 
described  in  the  flowchart  of  load  application  Fig.  4.3,  which  is  dis¬ 
cussed  later.  Once  the  solution  has  been  found  for  the  last  subincrement 
of  the  current  load  step,  the  overall  change  in  strain  for  the  wall 
segment  is  checked  to  see  if  this  change  is  too  large  for  the  number  of 
subincrements.  If  this  is  true,  the  solution  is  rejected  and  the  process 
returns  to  the  start  of  the  load  application,  with  the  values  of  the 
material  parameters  for  each  layer  recalled  from  the  end  of  the  last 
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load  step.  The  number  of  subincrements  is  increased  and  the  solution 
process  begins  again. 

The  method  of  solution  for  a  load  subincrement  is  described  in  the 
flowchart  of  load  application  shown  in  Fig.  4.3.  First  the  segment 
stiffness  matrix  is  compiled  from  the  results  of  the  last  trial  solution. 
The  program  loops  over  all  the  layers  calling  for  the  appropriate  con¬ 
stitutive  relations,  finds  the  stiffness  of  each  layer  and  integrates 
these.  Then  with  the  current  estimate  of  the  segment  strains,  the 
biaxial  strains  for  each  layer  are  found,  the  stresses  in  each  layer 
computed,  and  the  load  resistance  of  the  wall  segment  calculated.  The 
load  resistance  is  then  compared  to  the  applied  load.  If  the  difference 
is  small,  the  solution  has  been  found.  If  the  difference  is  large,  a 
correction  is  made  to  the  segment  strains,  and  the  loop  repeated. 

The  constitutive  relations  for  the  steel  layers  and  the  bilinear 
elastic  constitutive  relation  for  the  concrete  are  straight  forward. 
However,  the  nonlinear  concrete  constitutive  relation  is  more  complex 
and  is  described  in  Fig.  4.4.  This  relation,  if  specified,  is  called 
separately  for  each  concrete  layer  for  every  trial  solution.  The  solu¬ 
tion  for  the  layer  stiffness  and  stress  begins  by  again  checking  the  change 
in  real  strains.  If  this  change  is  too  large,  the  strain  change  is  sub¬ 
divided.  For  the  current  subdivision,  an  initial  estimate  is  made  of 
the  change  in  equivalent  uniaxial  strains,  Ae^  *,  using  the  values  of 
the  tangential  moduli,  E  ,  the  stresses,  o_^,  and  Poisson's  ratio,  u 
from  the  previous  subdivision.  This  initial  estimate,  Ac.  1  is  added 

1U 

to  the  previous  total  equivalent  uniaxial  strains,  e  ,  to  find  an 
estimate  of  the  new  total,  e.  Using  a,  and  c.  *,  the  failure  para- 
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meters  a.  1  and  z.  1  which  are  the  ultimate  stresses  for  the  current 

1C  1UC 

stress  ratio  and  the  equivalent  uniaxial  strains  at  ultimate  stresses 

for  the  current  stress  ratio,  are  calculated.  New  values  for  the  stress 

tangential  modulus  E_^ 1  and  Poisson’s  ratio  u  * ,  are  found  based  on 

z.  1 ,  a .  1  and  z .  1 .  A  second  estimate  of  Ac.  is  made,  based  on  E . 1 

1U  1C  iuc  iu  1 

and  u1.  This  new  estimate  Ac.  2  is  averaged  with  Ac.  1  to  find  Ac.  a. 

iu  iu  iu 

This  averaged  value,  Ae .  is  added  to  the  value  of  z.  from  the  last 

iu  iu 

subdivision  to  obtain  the  new  total  equivalent  uniaxial  strain,  E_^n. 
Based  on  a.*  and  z.  n,  the  failure  parameters  are  revised  to  a.  n  and 


iu 
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z .  n.  These  are  then  used  to  calculate  the  final  values  of  a.n,  E.n 
iuc  1  1 

and  un.  This  process  is  repeated  for  all  of  the  strain  subdivisions. 

In  total  then,  there  is  a  hierarchy  of  three  subdivision  loops. 

The  increment  generator  divides  each  load  into  load  steps.  This  number 
is  user  controlled.  The  second  subdivision  occurs  when  the  change  in 
strains  to  solve  a  load  step  is  too  large.  Here  the  load  step  is  sub¬ 
divided  into  load  increments.  The  final  increment  generator  is  within 
the  nonlinear  concrete  subroutine.  If  the  strain  change  for  the  solution 
of  a  load  increment  is  too  large,  then  this  strain  change  is  subdivided. 
If  however,  the  bilinear  concrete  equations  are  specified,  only  the 


first  subdivision  generator,  which  is  user  controlled,  is  operative 


4.3  Convergence  Difficulties 

Three  different  types  of  convergence  problems  may  be  encountered 
when  using  the  program.  The  most  common  type  of  convergence  difficulty 
involves  a  trial  solution  for  the  concrete  which  is  so  far  from  the 
correct  solution  that  subsequent  corrections  to  this  trial  lead  to 
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greater  error  in  the  solution.  The  end  result  of  this  incorrect 
solution  path  may  be  either  convergence  on  an  incorrect  solution  which 
is  in  equilibrium  with  the  loads,  or  no  convergence  at  all.  The 
three  main  levels  of  increment  generators  are  designed  to  minimize  the 
stray  movement  of  the  solution  path  by  keeping  the  strain  change  between 
solutions  small.  After  the  solution  of  each  load  increment,  the  strain 
change  is  calculated  and  if  this  change  is  too  large,  the  process 
returns  to  the  last  load  increment,  and  this  load  increment  is  sub¬ 
divided  into  load  subincrements . 

This  type  of  convergence  difficulty  was  most  often  encountered 
when  the  concrete  stress  is  moving  across  the  peak  of  the  stress-strain 
curve  or  in  other  words,  from  the  ascending  to  the  descending  branch. 
However  it  also  occurred  at  the  time  of  changes  in  the  stiffness  of  the 
reinforcing  bars  or  the  prestressing  strands.  This  type  of  convergence 
is  not  experienced  when  using  the  bilinear  elastic  concrete  relation. 

The  second  type  of  convergence  problem  which  may  be  encountered 
can  occur  with  either  the  bilinear  or  nonlinear  concrete  relations. 

This  problem  occurs  when  the  program  is  only  able  to  find  solutions 
which  bracket  the  correct  solution,  but  which  do  not  satisfy  the  accuracy 
requirements.  This  problem  occurs  most  frequently  when  a  layer  of  steel 
yields.  To  reduce  the  frequency  of  this,  and  also  the  first  type  of 
convergence  problem,  an  under-relaxation  factor  is  used  which  is  applied 
to  and  reduces  the  calculated  strain  corrections.  The  under-relaxation 
factor  is  calculated  to  maintain  the  strain  corrections  below  a  set  value. 
If,  however,  a  solution  which  meets  the  accuracy  requirements  of  0.1% 
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is  not  found  within  90  attempts,  the  accuracy  requirement  is  reduced 
to  5%,  and  if  this  limit  is  satisfied,  the  solution  is  accepted  and  the 
next  load  increment  is  attempted.  The  determination  of  the  accuracy  of 
a  solution  is  discussed  in  the  next  section. 

If  strains  rather  than  forces  are  specified,  none  of  these 
difficulties  will  occur.  If  a  segment  response  curve  contains  a  portion 
of  overall  negative  stiffness,  then  again  no  difficulties  will  occur 
if  strains  are  specified.  However  if  forces  are  specified,  then  portions 
of  negative  stiffness  for  the  wall  segment  can  cause  difficulties.  To 
force  the  segment  to  accept  increasing  loads,  after  20  strain  corrections 
have  been  applied,  the  strains  are  constrained  to  move  in  the  same 
direction  as  the  loads.  Thus,  if  the  axial  force  in  Direction  1  is 
increasing,  then  the  strains  in  this  direction  are  not  allowed  to  decrease. 
This  mechanism  can  be  successful  in  jumping  small  dips  in  the  response 
curve . 

A  fourth  convergence  problem  which  always  occurred  after  the  con¬ 
crete  response  had  moved  into  the  post  ultimate  region  of  the  stress- 
strain  curves,  involved  the  degeneration  of  the  stress-strain  curve  in 
the  minor  direction.  This  problem  has  been  effectively  eliminated.  The 
technique  employed  to  overcome  this  problem  is  dealt  with  in  Section 


3.3.7. 


4.4  Solution  Accuracy 

The  accuracy  of  a  solution  is  approximated  by  the  following  equation 


(4.4.1) 
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where  AQ_^  is  the  unbalanced  load  at  degree  of  freedom  i,  is  the 
specified  load  at  degree  of  freedom  i,  m  is  the  number  of  degrees  of 
freedom  where  a  force  rather  than  a  deformation  is  specified,  and  the 
result  A  is  the  calculated  accuracy  of  the  solution.  Only  those  degrees 
of  freedom  at  which  forces  are  specified  are  included  in  this  equation 
so  that  (AQ^  /  F_^)  is  set  to  zero  for  those  degrees  of  freedom  at  which 
deformations  are  specified.  The  accuracy  limit  which  A  must  satisfy 
for  solutions  determined  in  ninety  trials  or  less  is  .001.  For  solutions 
found  between  the  ninety-first  and  the  one  hundred  and  twentieth  trial, 

A  must  be  less  than  or  equal  to  .05. 

4.5  Layering  Technique 

The  concrete  wall  segment,  when  subjected  to  a  curvature,  does  not 
react  uniformly  across  the  wall  thickness  because  of  the  imposed  strain 
gradient.  In  addition,  the  presence  of  the  steel  layers  causes  further 
differences  in  the  response  of  the  segment,  from  that  of  a  plain  concrete 
section.  To  account  for  these  effects,  the  wall  segment  has  been  sub¬ 
divided  into  layers  through  the  wall  segment  as  shown  in  Fig.  4.6. 

Within  each  layer,  the  material  is  assumed  to  have  uniform  properties 
and  individual  layers  are  not  influenced  by  adjoining  layers.  By  then 
integrating  across  all  of  the  concrete  and  steel  layers,  the 
composite  response  of  the  wall  segment,  including  any  nonlinear  results, 
will  be  assessed.  The  values  of  the  elements  in  the  matrix  [c]  of  Eq . 
2.3.25  and  in  the  vector  of  Eq.  2.3.19  are  found  by  calculating  separately 
for  each  concrete  and  steel  layer,  the  values  of  a^,  Z  and  3f  /  3e^. 

The  integration  of  these  values  through  the  entire  thickness  of  the  wall 
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gives  the  complete  wall  properties. 

This  integration  is  organized  into  two  parts.  First,  the  wall 
segment  is  modelled  ignoring  the  steel.  Second,  the  steel  is  added. 

The  integration  of  the  plain  concrete  follows  the  stated  technique 
without  complication.  The  advantage  of  the  technique  is  that  the  number 
of  the  plain  concrete  layers,  and  their  thicknesses  are  uniform,  not 
being  affected  by  the  steel.  The  steel  layers,  at  any  location  in  the 
segment,  and  of  any  thickness,  are  simply  added  into  the  sum.  To  account 
for  the  concrete  which  was  displaced  by  the  steel,  but  has  been  included 
in  the  integrated  properties,  a  negative  concrete  layer  is  associated  with 
each  steel  layer.  The  results  obtained  from  this  negative  concrete  layer 
are  subtracted  from  the  integration  of  the  wall  segment.  Because  of  the 
unidirectional  contribution  of  the  steel,  the  negative  concrete  layer 
subtracts  only  in  the  same  direction  as  the  steel.  The  materials 
quantities  of  the  negative  concrete  layer,  however,  are  calculated  in  a 
manner  consistent  with  the  plain  concrete  layers. 

4.6  Example  Problem 

The  example  segment  which  shall  be  examined  is  taken  from  the 
elastic  stress  analysis  report  by  Epstein  and  Murray  (1976) ,  henceforth 
referred  to  as  the  EM  report,  and  the  segment  is  labelled  in  this  report 
as  UD1.  This  doubly-curved  segment  is  located  in  the  upper  dome,  next 
to  the  ring  beam  of  the  Gentilly-type  structure  which  in  part  consists 
of  the  outside  dome  and  an  inner  dome  which  together  form  a  water 
reservoir.  Although  in  the  actual  structure,  the  dome  is  thickened 
where  it  meets  the  ring  beam,  in  this  analysis  it  is  assumed  that  this 
thickening  is  not  present  and  the  analysis  will  ignore  any  extra  strength 
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that  the  thickening  may  provide.  It  does  not,  therefore,  accurately 
reflect  the  response  which  would  be  expected  in  the  prototype  structure. 

4.6.1  Segment  Description 

The  geometry  of  the  segment  and  the  locations  of  the  reinforcing 
steel  and  prestressing  steel  are  shown  in  Fig.  4.7.  The  area  of  the  steel 
layers  and  the  prestressing  forces  are  listed  in  Table  4.1.  The  rein¬ 
forcing  layers  have  a  yield  stress  of  60  ksi  and  have  a  Young’s  modulus 
of  29000  ksi.  The  prestressing  steel  has  an  ultimate  strength  of  255  ksi 
with  an  initial  modulus  of  29600  ksi. 

The  concrete  strength  in  uniaxial  compression  for  this  example 
was  assumed  to  be  5  ksi.  Using  this  value  and  the  A.C.I.  recommendations 
(American  Concrete  Institute,  1971)  the  initial  concrete  modulus  was 
calculated  to  be  4030  ksi.  The  tensile  strength  of  0.424  ksi  was 
computed  using  6  /  f '  .  An  initial  Poisson's  ratio  was  selected  to  be 
0.15.  The  default  value  provided  in  the  program  was  used  for  e  .  The 

value  of  e.  was  input  as  0.000126  or  1.2  f'  /  E  .  For  both  the  con- 
lot  to 

Crete  and  steel,  the  coefficient  of  thermal  expansion  was  assumed  to  be 
0.65  x  10  in.  /  in.  F  . 

4.6.2  Loading  Description 

To  determine  the  forces  and  moments  to  which  the  segment  will  be 
subjected,  the  containment  vessel  was  analysed  using  a  linear  elastic 
finite  difference  program  based  on  the  classical  thin  shell  equations 
(Epstein  and  Murray,  1976)  .  From  this  program  the  forces  and  moments 
resulting  from  dead  loads,  prestressing  loads,  and  live  loads  were  found. 
Table  4.2  lists  the  total  applied  load  vectors  to  which  the  dome  segment 
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was  subjected  for  each  load  step.  As  part  of  modelling  the  construction 
sequence,  a  small  dead  load  was  applied  before  prestressing.  After 
prestressing,  a  period  of  time  will  pass  before  the  live  loads  due  to 
overpressure  are  applied,  and  therefore  a  creep  step  is  provided.  The 
creep  factor  is  set  to  2.50.  The  thermal  loads  are  simulated  by  a 
linear  temperature  gradient  varying  from  +26°F  on  the  inner  surface  to 
-50°F  on  the  outer  surface,  the  same  temperature  distribution  as  used 
in  the  EM  report.  For  this  example,  it  was  assumed  that  the  segment 
would  be  restrained  against  curvatures  resulting  from  the  temperature 
gradient,  but  that  axial  movement  would  not  be  restricted.  In  other 
words,  the  temperature  change  would  affect  the  in-plane  axial  deforma¬ 
tions  and  out-of-plane  moments,  but  would  not  change  the  axial  forces 
or  curvatures.  The  live  loads  are  the  last  to  be  imposed.  The  moments 
and  axial  forces  resulting  from  a  unit  of  internal  pressure  are  listed 
in  the  EM  report  and  were  calculated  using  the  classical  shell  analysis. 
It  was  then  assumed  that  the  ratio  of  axial  force  to  curvature  resulting 
from  the  unit  of  internal  pressure  would  remain  constant  up  to  the 
ultimate  load  of  the  segment.  By  specifying  the  live  load  to  be  of  this 
type,  as  was  used  in  the  EM  report,  the  ratio  of  Nj  /  Mj  and  N2  /  M2 
can  be  input  to  COWSAC  and  the  necessary  curvatures  are  calculated  using 
a  linear  elastic  uncracked  section. 

4.6.3  Discussion  of  Results 

Figure  4.8  contains  the  axial  load-axial  strain  curves  for  both 
directions  one  and  two.  The  various  points  labelled  1.1  -  1.5  and 
2.1  -  2.5  are  the  results  from  the  various  applied  loads  up  to  and 
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including  the  live  loads.  Points  1.1  and  2.1  represent  the  segment 
response  to  the  dead  load  in  directions  1  and  2  respectively.  Points 
1.2  and  2.2  describe  the  response  to  the  prestressing.  The  prestressing 
forces  are  internal  to  the  segment  but  the  external  axial  forces  are  not 
constant  therefore  both  the  strains  e  ^mS  and  C2mS  an<^  the  applied  forces 
Nj  and  N2  change.  Points  1.3  and  2.3  show  the  results  of  the  creep  step. 
During  this  operation,  the  loads  are  assumed  to  remain  constant  and 
therefore  only  the  strains  change.  Locations  1.4  and  2.4  exhibit  the 
results  of  the  temperature  changes.  Since  the  axial  degrees  of  freedom 
are  not  restrained  against  thermal  movement,  the  axial  forces  remain 
constant  while  the  axial  strains  change. 

After  the  solution  to  the  thermal  effects  has  been  determined, 
the  live  loads  are  applied.  As  previously  stated,  the  linear  elastic 
response  of  containment  vessel  to  a  unit  of  internal  pressure  was  computed 
by  the  authors  of  the  EM  report,  who  then  assumed  that  the  axial  forces 
and  the  curvatures  which  resulted  from  this  analysis  would  vary  in  a 
linear  manner  to  higher  internal  pressures.  This  loading  path  was  also 
employed  here  to  obtain  the  results  shown  for  this  example. 

The  solid  lines  of  Fig.  4.8  trace  the  load-strain  curves  for  the 
segment  in  both  directions,  for  a  set  of  loads  consistent  with  the  above 

ms 

hypothesis.  The  direction  labelled  N 1  -  c ^  starts  at  an  external  load 
level  near  zero  and  then  moves  into  the  tension  zone .  The  response  is 
linear  up  to  approximately  200  kips  whereupon  the  segment  begins  to 
soften  gradually.  If  the  segment  were  loaded  uniformly  with  no  moment, 
then  it  would  be  expected  that  this  initial  linear  behavior  would  persist 
up  to  the  cracking  load,  or  about  450  kips.  The  overall  segment  stiffness 
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would  then  change  abruptly  after  the  concrete  had  cracked  in  tension. 

This  type  of  behavior  is  demonstrated  by  Segment  UD3  in  the  EM  report. 
However,  the  presence  of  the  large  moment  in  direction  1  of  this  example 
segment,  UD1,  causes  some  of  the  concrete  to  crack  at  lower  loads  and 
the  sudden  break  in  the  load-strain  curve  that  would  be  expected  under 
these  conditions  is  replaced  with  the  smooth  transition  shown  in  Fig. 

4.  8. 

ms 

In  the  second  direction,  labelled  N2  -  £2  »  the  segment  begins  in 

tension  but  then  moves  into  compression.  A  large  moment  is  present  in 
this  direction  at  lower  loads,  but  the  applied  moment  disappears  as 
the  compressive  axial  force  increases.  Since  the  applied  axial  loads 
are  well  under  the  capacity  of  the  segment,  the  response  curve  in  this 
direction  is  linear. 

Also  shown  in  Fig.  4.8  are  the  response  curves  in  both  directions 
that  are  obtained  using  the  same  loading  conditions  but  with  the  bilinear 
elastic  model.  These  results  are  designated  with  the  dashed-dotted  lines, 
and  are  the  results  which  would  be  expected  from  the  EM  report.  These 
curves  indicate  that  there  are  no  significant  differences  between  the 
results  obtained  from  the  nonlinear  or  the  bilinear  concrete  constitutive 
equations  as  described.  However,  the  results  obtained  from  this  particular 
problem  are  not  general  in  that  the  very  high  moment  reduces  the  signifi¬ 
cance  of  the  tensile  behavior  of  the  concrete.  The  uniqueness  of  these 
results  is  demonstrated  by  the  results  in  Chapter  5,  which  do  not  exhibit 
a  close  comparison  between  the  nonlinear  and  bilinear  solutions  for  the 
test  segments.  In  fact,  the  EM  report  showed  no  differences  in  response 
for  UD1  when  the  tensile  strength  of  the  concrete  was  reduced  to  zero 
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from  6  /  f '  . 

c 

The  moment- curvature  diagram  for  Segment  UD1  when  subjected  to 
the  above  loading  sequence  is  illustrated  in  Fig.  4.9.  The  points 
labelled  1.1  to  1.5  and  2.1  to  2.5  are  similar  to  those  in  the  axial 
force-strain  diagram.  Points  1.1  and  2.1  show  the  moment  and  curvature 
in  directions  1  and  2  respectively,  after  the  application  of  dead  load. 
With  the  application  of  the  prestressing  forces,  a  large  moment  is  placed 
on  direction  1  and  a  small  moment  is  placed  on  direction  2.  These 
moments  are  not  directly  the  result  of  the  prestressing  but  are  applied 
moments  which  the  finite  difference  analysis  of  the  EM  report  predicts 
will  accompany  the  prestressing.  The  results  of  the  prestressing  is 
indicated  by  the  points  labelled  1.2  and  2.2.  After  the  prestressing, 
the  segment  is  allowed  to  creep.  This  step  is  shown  by  locations  1.3 
and  2.3.  The  applied  moments  and  M2  are  assumed  to  remain  constant, 
but  due  to  creep,  (j>i  and  <£2  change  from  the  previous  step.  During  the 
calculation  of  the  thermal  effects  however,  it  was  assumed  that  the  seg¬ 
ment  was  restrained  so  that  the  curvatures  were  kept  constant  and  the 
moments  allowed  to  vary.  The  results  of  the  linear  temperature  variation 
are  shown  at  points  1.4  and  2.4.  The  effect  of  increasing  the  temperature 
on  the  inner  surface  and  decreasing  the  temperature  on  the  outer  surface 
was  a  large  negative  moment  in  both  directions . 

The  locations  where  live  loading  begins  are  indicated  by  points 
1.5  and  2.5  which  differ  from  the  previous  solution  at  1.4  and  2.4.  This 
discrepancy  is  due  to  the  fact  that  some  of  the  outside  concrete  layers 
unloaded  after  they  reached  their  maximum  tensile  stress  and  were  on  the 
declining  branch  of  the  stress-strain  curve.  This  unloading,  the  defini- 
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tion  of  which  is  given  in  Chapter  3,  is  the  consequence  of  the  live 
load  moment  opposing  the  moment  after  the  application  of  the  thermal 
effects.  When  unloading  is  detected  in  a  concrete  layer  which  is  on 
the  declining  branch  of  a  tensile  stress-strain  curve,  the  stresses  and 
stiffness  in  that  concrete  layer  drop  to  zero  and  remain  there  until 
that  layer  goes  into  compression.  These  layers  can  no  longer  resist  any 
tensile  loads.  Since  these  layers  now  do  not  contribute  strength  to  the 
segment,  a  sudden  reduction  in  moment  resistance  is  experienced. 

Beginning  at  locations  1.5  and  2.5,  the  solid  lines  trace  the 
response  of  the  segment  to  the  live  loads.  The  curve  of  -  <£1  begins 
with  an  initial  section  of  increasing  stiffness.  The  stiffness  continues 
to  increase  until  all  of  the  cracked  layers  are  closed  and  are  accepting 
stresses.  From  this  point,  the  curve  is  linear  up  to  the  region  where 
the  concrete  layers  behave  inelastically .  Here  the  slope  of  the  curve 
declines  to  zero,  at  which  point  the  maximum  moment  is  reached.  Since 
the  loading  is  curvature  controlled,  the  moment  now  begins  to  decrease 
with  increasing  curvature. 

The  behavior  of  the  moment-curvature  trace  in  the  second  direction 
is  similar  to  that  in  the  first  direction.  Again,  there  is  the  jump  from 
the  thermal  load  solution  to  the  start  of  live  loading,  dictated  by  the 
unloading  model.  After  this,  the  moments  increase  linearly  up  to  the 
point  where  loading  stops,  at  a  value  comparable  to  the  maximum  moment 
in  the  first  direction.  Since  the  two  directions  of  Segment  UD1  are 
similar  in  details  of  construction  any  apparent  difference  in  moment 
capacity  between  directions  1  and  2  is  a  result  of  the  axial  loading. 

Also  shown  in  Fig.  4.9  are  the  response  curves  in  both  directions 
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from  the  bilinear  elastic  analysis  using  the  same  loading  conditions. 

These  results  are  denoted  by  the  dashed-dotted  lines.  Again,  as  for  the 
axial  force-strain  diagram,  the  bilinear  and  nonlinear  solutions  compare 
well.  The  only  characteristic  not  displayed  by  the  linear  analysis  is 
the  jump  from  the  thermal  solution  to  the  start  of  live  loading.  However 
this  jump  could  not  be  expected  from  the  bilinear  solution  as  this  con¬ 
crete  model  has  no  post  ultimate  portion,  and  therefore  no  unloading 
from  this  region.  However  in  direction  1,  the  nonlinear  solution 
exhibits  extra  moment  capacity  after  a  curvature  of  approximately 
.00002  in. /in. 

Some  of  the  characteristics  displayed  by  the  M  —  <f>  plots  can  be 
further  understood  by  examining  the  stress  path  of  the  concrete. 

Figure  4.10  shows  the  concrete  stress  path  for  the  concrete  on  the 
inner  surface  of  the  segment,  and  Figure  4.11  for  the  concrete  on  the  outer 
surfaces.  Locations  1-4  indicate  the  solutions  to  the  dead,  prestressing 
creep  and  thermal  loads.  The  solid  line  from  positions  4  to  5  marks  the 
stress  path  followed  during  the  application  of  the  live  load.  For  the 
concrete  on  the  inner  surface  of  Segment  UD1,  Fig.  4.10,  the  direction 
of  this  path  results  in  an  initially  constant  a 2,  as  o\  moves  toward  the 
a 2  axis.  After  01  becomes  positive,  the  01-02  stress  path  touches  the 
failure  envelop  and  01  decreases  to  its  minimum  value. 

The  consequence  of  unloading  on  the  declining  branch  of  a  stress- 
strain  curve  which  is  discussed  in  Chapter  3,  are  shown  in  Fig.  4.11. 

After  the  creep  solution  has  been  found,  shown  at  location  3  in  Fig.  4.11, 
the  outside  layer  of  concrete  has  positive  stiffness.  When  the  thermal 
strains  are  added,  the  stresses  hit  the  failure  envelope  and  start  to 
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decline.  Both  directions  are  now  on  the  declining  branch  of  their 
a  -  curves.  Under  the  application  of  internal  pressure  the  applied 

moment  now  reverses  and  the  outer  concrete  layers  begin  to  unload.  The 
response  of  the  nonlinear  concrete  constitutive  model  to  unloading  from 
the  declining  branch  is  to  set  the  concrete  stresses  to  zero,  and  so  the 
next  point  in  Fig.  4.11  after  point  4  is  at  the  origin  of  the  stress 
plot.  As  the  loading  continues,  the  real  strains  continue  to  decrease 
and  eventually  the  cracks  close  in  direction  1,  and  this  direction 
picks  up  compressive  stresses,  and  the  stress  path  begins  to  move  along 
the  cji  axis.  Further  in  the  loading,  the  cracks  transverse  to  direction 
2  close  and  this  direction  accepts  compressive  stresses.  At  this  point 
the  stress  path  leaves  the  oi  axis.  As  the  tensile  axial  load  increases, 

the  stress  a^,  again  moves  toward  the  tensile  region.  However,  because 
of  the  prior  cracking,  this  layer  of  concrete  cannot  accept  any  tensile 
stresses  and  the  stress  path  is  restricted  to  the  0£  axis. 

Figure  4.12  details  the  depth  of  open  cracks  transverse  to  direction 
1  as  the  axial  live  load  Nj  increases.  At  the  start  of  live  loading,  the 
open  cracks  on  the  exterior  of  the  segment  exist  through  12  %  of  the  seg¬ 
ment  thickness.  As  the  live  load  increases,  these  cracks  close  until  the 
live  axial  load  reaches  75  kips,  at  which  point  no  open  cracks  exist  on 
the  outside  of  the  segment.  At  a  line  axial  load  of  200  kips,  cracks 
begin  to  appear  on  the  interior  side  of  the  segment.  The  depth  of  these 
open  cracks  increases  with  increasing  live  load  until  this  load  reaches 
575  kips,  at  which  point  open  cracks  penetrate  through  the  wall  segment. 

Also  shown  in  Fig.  4.12  are  the  cracking  results  obtained  by  Epstein  and 
Murray  for  Segment  UD1 .  The  results  by  the  two  different  analysis  agree, 


the  difference  usually  being  less  than  50  kips. 
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:  e 1  Layer 

Area  (in2) 

An 

1.1 

A12 

1.1 

A21 

1.27 

A22 

1.69 

?! 

2.33 

P2 

2.33 

TABLE  4.1 

Steel  Layer 

Prestress  Force  (kips) 


374.0 

200.0 

Details  -  Segment  UD1 
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MAIN  OF  COWSAC 


Fig.  4.1  Flowchart  of  Main  of  COWSAC 
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SUBROUTINE  WORK 


Fig.  4.2  Flowchart  of  Load  Control 


76 


SUBROUTINE  WORK  1 


*  For  bilinear  concrete  material,  call  SUBROUTINE  LINELC 

For  nonlinear  concrete  material,  call  SUBROUTINE  SAENZ 


Fig.  4.3  Flowchart  of  Load  Application  and  Solution 
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SUBROUTINE  SAENZ 


Fig.  4.4  Flowchart  of  Nonlinear  Concrete  Constitutive  Relation 
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SUBROUTINE  LINELC 


Fig.  4.5 


Flowchart  of  Bilinear  Concrete  Constitutive  Relation 
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Fig.  4.6  Wall  Segment  Showing  Orientation  of  Layers 
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Fig.  4.7  Segment  UDl 
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Fig.  4.8  Applied  Axial  Load  vs.  Strain  -  Segment  UD1 
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Result  Co-ordinate 
Position 

Dead  Load 
Prestressing 
Creep 

Temperature 
Live  Loads 


-4  -3  -2  -1  0 

Stress  in  Direction  1  (ksi) 


Fig.  4.10  Concrete  Biaxial  Stress  Path  on  Inner  Surface  of  Segment  UD1 


Stress  in  Direction  2  (ksi) 
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Fig.  4.11  Concrete  Biaxial  Stress  Path  on  Outer  Surface  of 


Segment  UD1 
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5.  APPLICATIONS 


5.1  Comparison  of  Program  with  Wall  Segment  Tests 

The  ultimate  test  of  the  accuracy  and  adequacy  of  a  computer  pro¬ 
gram  and  the  theory  on  which  it  is  based  is  the  comparison  between  test 
results  and  predicted  results.  To  test  this  computer  model,  the  program 
output  was  compared  with  the  observations  from  eight  wall  segment  tests, 
one  including  imposed  moments  (Rizkalla,  et  al,  1979) .  The  wall  segments, 
of  common  width  and  length,  were  cast  in  two  different  thicknesses,  with 
a  varying  amount  of  reinforcing  steel,  and  in  some  specimens,  with  pre¬ 
stressing  steel.  Three  different  load  ratios  were  used.  The  two  cross- 
sections  of  interest  for  Specimen  No.  1  are  shown  in  Fig.  5.1  and  5.2. 
Tables  5.1  and  5.2  contain  a  summary  of  the  pertinent  information  of  all 
the  wall  segments.  Of  the  data  presented  in  this  table,  it  should  be 
noted  that  the  value  of  f’  for  each  specimen  was  obtained  from  a 
Brazilian  tensile  test.  The  best  value  to  be  used  as  the  ultimate  ten¬ 
sile  strength  for  the  test  segments  was  found  to  be  0.6  f'  (Murray,  et  al, 
1978).  This  value  shall  be  accepted  in  the  following  analyses.  The 
value  for  the  initial  elastic  modulus  of  the  concrete  was  assumed 
(Rizkalla,  et  al,  1979)  to  be  3800  ksi  for  all  the  specimens.  Similarly, 
the  initial  value  of  Poisson's  ratio  was  fixed  at  0.2. 

The  tests  which  were  conducted  on  these  specimens  were  primarily 
tension  tests,  since  tensile  live  loads  were  placed  on  the  wall  segments 
and  all  of  the  walls  failed  in  tension.  Although  this  provides  sufficient 
data  to  verify  the  degrading  stiffness  used  for  concrete  in  tension 
data  are  not  available  to  check  the  failure  curves  of  the  concrete  in 
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tension-compression  and  compression-compression,  or  to  accurately  assess 
the  declining  branches  of  the  Saenz  stress-strain  curve  in  compression. 
The  presence  of  the  prestressing,  however,  allows  the  opportunity  to 
check  the  program  in  the  compression-compression  and  tension-compression 
quadrants  at  levels  of  stress  below  the  maximum,  since  after  the  applica¬ 
tion  of  the  prestress,  the  wall  segment  is  subjected  to  biaxial  compres¬ 
sion.  As  the  tensile  loads  increase,  the  stresses  move  into  the  tension- 
compression  quadrant,  and  finally  into  the  biaxial  tension  quadrant. 

This  path  is  typified  by  the  stress  path  of  Specimen  No.  1  shown  in 
Fig.  5.3. 

To  produce  the  tensile  loads  on  the  segment,  the  specimen  was 
attached  to  the  loading  machine  by  fixing  the  loading  heads  onto  the 
ends  of  the  prestressing  strands  and  reinforcing  bars.  There  was  no 
direct  connection  between  the  testing  apparatus  and  the  concrete.  All 
forces  had  to  be  transferred  from  the  loading  jacks  to  the  concrete 
through  the  steel  elements.  Due  to  the  assumptions  made  in  the  analysis, 
this  system  of  loading  cannot  be  modelled  and  it  will  be  assumed  that  the 
entire  face  of  the  segment,  rather  than  only  the  steel  layers,  will  be 
loaded . 

In  these  tests,  the  major  contribution  to  strength  is  made  by 
the  reinforcing  and,  if  present,  the  prestressing  steel.  For  this 
reason,  the  specimen  axial  force-strain  response  plots  contain  curves 
which  represent  the  capacity  at  a  given  strain,  of  the  steel  only. 
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5.1.1  Wall  Segments  with  Prestressing 

Segments  Nos.  1,  3,  8  and  12  contained  prestressing  in  both  direc¬ 
tions  and  were  tested  at  load  ratios  of  2:1,  1:1,  2:1  and  2:1  respectively, 
the  latter  specimen  including  a  moment  in  one  direction.  Segments  Nos. 

5  and  6  contained  prestressing  in  one  direction  only  and  were  both  tested 
at  load  ratios  of  1:0. 

5. 1.1.1  Wall  Segment  No.  1 

This  segment  was  tested  at  a  load  ratio  of  2:1  with  the  larger 
imposed  force  being  in  direction  1,  as  designated  in  Table  5.1,  which 
is  the  direction  containing  the  greater  area  of  prestressing  tendons. 

The  tensile  strength  of  the  concrete  was  found  by  the  Brazilian  test 
to  be  490  psi.  The  effective  tensile  strength  was  taken  to  be  six-tenths 
of  this  value  or  294  psi. 

Up  to  the  initial  cracking  load  of  the  concrete  in  either  direction, 
the  program  results  predict  a  stiffer  response  than  observed,  as  shown 
in  Fig.  5.4.  This  however,  is  highly  dependent  upon  the  chosen  value 
for  the  initial  elastic  modulus  of  the  concrete.  More  accurate  data  for 
this  variable  may  provide  closer  agreement.  In  the  major  direction 
(direction  1),  the  model  predicts  a  higher  cracking  load  than  that 
observed,  and  that  additional  concrete  strength  continues  to  influence 
the  response  in  the  major  direction  until  the  concrete  strength  is 
exhausted.  However,  the  stiffness  of  the  model  is  close  to  that  of  the 
test  specimen  after  cracking.  In  the  minor  direction,  the  predicted 
curve  is  very  close  to  the  observed  points. 
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5. 1.1. 2  Wall  Segment  No.  3 

Specimen  No.  3,  containing  the  same  steel  and  prestressing  as 
Specimen  No.  1,  was  tested  at  a  load  ratio  of  1:1  up  to  375  kips,  after 
which  direction  1  was  loaded  further,  keeping  the  load  on  direction  2 
constant.  The  tensile  strength  of  the  concrete  was  426  psi,  giving  an 
effective  strength  of  256  psi. 

Up  to  cracking  there  is  good  agreement  between  the  response  curves 
of  the  model  and  the  test  specimen,  which  are  shown  in  Fig.  5.5.  In 
direction  2,  a  close  correlation  exists  between  the  predicted  and  observed 
curves,  the  model  predicting  an  extra  load  capacity  of  about  20  kips  at 
a  given  strain.  Again  however,  the  specimen  stiffness  of  the  model  after 
cracking  is  very  close  to  the  observed  results.  In  direction  1,  there 
is  a  slow  change  in  stiffness  at  the  ’cracking  load',  and  this  is  pre¬ 
dicted  by  the  model.  After  the  change  in  loading  ratio,  the  response 
curve  of  the  model  is  more  flexible  than  that  shown  by  the  test  specimen. 

5.1.1. 3  Wall  Segment  No.  8 

This  specimen  again  contains  the  same  steel  and  prestressing 
arrangement  as  Specimens  Nos.  1  and  3,  and  was  tested  at  a  load  ratio 
of  2:1,  as  was  Segment  No.  1.  The  tensile  strength  of  the  concrete 
was  424  psi  resulting  in  an  effective  tensile  strength  of  254  psi. 

Up  to  cracking,  the  agreement  between  the  model  and  the  specimen 
is  good  in  both  directions,  the  choice  of  the  initial  elastic  modulus  of 
the  concrete  being  the  most  significant  factor.  In  direction  1,  the 
cracking  load  of  the  model  and  test  specimen  are  similar,  with  the 
response  of  the  model  about  20  kips  below  the  curve  of  the  specimen. 
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These  curves  are  shown  in  Fig.  5.6.  Both  curves  have  similar  slopes 
after  cracking.  In  direction  2,  the  agreement  of  the  cracking  load  and 
of  the  stiffness  after  cracking  between  the  model  and  the  test  specimen 
is  very  good. 

5. 1.1. 4  Wall  Segment  No.  5 

Specimen  No.  5  contains  prestressing  steel  only  in  one  direction  and 
it  was  loaded  only  in  the  prestressed  direction,  for  a  load  ratio  of 
1:0.  The  concrete  tensile  strength  was  420  psi,  the  effective  tensile 
strength,  252  psi. 

The  model  and  the  test  again  provide  reasonable  agreement,  with 
the  model  predicting  load-carrying  capabilities  of  about  20  kips  less 
after  cracking.  The  stiffness  after  cracking  of  the  model  and  the  test 
segment  are  similar.  The  results  from  this  segment  are  shown  in  Fig.  5.7. 

5. 1.1. 5  Wall  Segment  No.  6 

Segment  No.  6,  like  Specimen  No.  5,  contains  prestressing  in  only 
one  direction,  but  more  reinforcing  steel  is  present.  The  load  is  again 
applied  only  in  the  direction  of  the  prestressing.  The  measured  tensile 
strength  was  361  psi  for  an  effective  strength  of  217  psi. 

The  model  and  the  test  segment  agree  reasonably  well,  both  sharing 
a  common  total  stiffness  after  cracking,  but  with  the  model  predicting 
a  curve  about  30  kips  above  that  measured  (See  Fig.  5.8). 

5. 1.1.6  Wall  Segment  No.  12 

This  wall  segment,  similar  in  construction  to  Segments  Nos.  1,  3 
and  9,  has  axial  loads  applied  in  the  ratio  of  2:1  and  additionally  a 
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moment  across  the  wall  face  orthogonal  to  the  minor  direction.  The 
applied  load  on  this  specimen  is  complicated  by  the  method  in  which 
forces  are  transferred  from  the  testing  machine  to  the  wall  segment. 

Since  this  transfer  occurs  through  the  prestressing  steel  and  the  rein¬ 
forcing  steel,  strain  measurements  on  the  exposed  lengths  of  the  steel 
were  used  to  compute  the  actual  loading  history  in  the  direction  with 
moment.  The  loading  sequence  for  this  analysis  is  listed  in  Table  5.3, 
and  it  should  be  noted  that  the  curvature  in  the  first  direction,  <f>  i , 
was  assumed  to  be  constant  at  zero.  For  the  purposes  of  discussion, 

Face  A  shall  be  the  inside  face  of  the  segment,  so  that  the  applied 
positive  moment  produces  tensile  strains  on  this  face. 

The  presentation  of  the  test  observations  and  the  results  of 
both  the  linear  and  nonlinear  analysis  are  shown  in  Fig.  5.9  and  5.10. 

On  Face  A,  the  predicted  response  curve  of  the  nonlinear  analysis  in 
direction  1  lies  above  that  observed  during  the  test.  It  appears  from 
the  test  results  that  the  concrete  is  not  contributing  to  the  segment 
stiffness  in  this  direction,  as  the  test  observations  follow  the  response 
of  the  linear  elastic  tension  cut-off  curve  after  cracking.  In  direction 
2,  the  curve  produced  by  the  nonlinear  analysis  closely  follows  the  test 
observations  up  to  250  kips.  After  this  point,  at  a  strain  of  approximate¬ 
ly  .002  in. /in.,  the  steel  in  the  analytical  model  yields  and  the  stiff¬ 
ness  of  the  section  is  reduced.  This  reduction  in  stiffness  did  not 
occur  in  the  test  specimen. 

On  Face  B  of  Specimen  No.  12,  the  predicted  load-strain  response  of 
the  nonlinear  analysis  is  close  to  the  response  of  the  test  specimen  for 
direction  1.  However,  as  in  the  curves  of  Face  A,  direction  2,  the 
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analytical  model  softens  at  the  yield  strain  of  .002  while  the  test  spe¬ 
cimen  does  not.  There  is  close  agreement  again  in  direction  2,  up  to  a 
load  of  270  kips.  Beyond  this  point,  the  predicted  response  diverges 
from  the  actual  response. 

The  close  agreement  between  the  analytical  results  and  the  test 
results  in  direction  2  for  both  Face  A  and  Face  B  indicates  that  the 
predicted  moment-curvature  relationship,  M2“<j>2>  should  correspond 
accurately  with  the  observed  results.  Since  zero  curvature  in  direction 
1  was  imposed  in  the  loading  conditions,  the  results  of  the  nonlinear 
analysis  show  the  response  curves  in  direction  1  on  Face  A  and  Face  B 
to  be  the  same.  However,  the  test  results  indicate  the  presence  of  a 
strain  gradient  in  this  direction.  Reanalysis  of  the  section  setting 
Mi  equal  to  zero,  as  opposed  to  setting  <f>i  equal  to  zero  during  loading, 
does  not  significantly  alter  the  results.  The  strain  gradient  from 
Face  A  to  Face  B  in  direction  1  must  then  have  been  produced  by  some 
detail  of  the  testing  procedure. 

5.1.2  Wall  Segments  Without  Prestressing 
5.1.2. 1  Wall  Segment  No.  4 

Specimen  No.  4  was  tested  at  a  load  ratio  of  1:1,  and  the  computed 
effective  strength  of  the  concrete  was  322  psi.  The  comparison  between 
the  model  and  the  specimen  is  poor.  The  cracking  load  was  measured  to 
be  about  50  kips  whereas  the  model  predicts  a  cracking  load  of  100  kips. 
This  difference  in  cracking  loads  of  50  kips  is  equivalent  to  approxi¬ 
mately  150  psi  in  concrete  stress.  After  cracking,  the  concrete  in  the 
test  specimen  maintains  its  strength  in  one  direction,  as  the  observed 
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load-strain  curve  remains  parallel  to  the  steel  only  curve,  while  the 
strength  decreases  in  the  orthogonal  direction.  Both  curves  eventually 
cross  under  the  steel  only  curve.  These  curves  are  shown  in  Fig.  5.11. 

5. 1.2. 2  Wall  Segment  No.  7 

Specimen  No.  7  was  tested  at  a  load  ratio  of  1:1.  The  effective 
concrete  strength  was  calculated  to  be  203  psi.  In  contrast  to  the  pre¬ 
vious  segments,  this  specimen  was  constructed  to  a  thickness  of  15.75 
inches . 

This  segment  also  compares  poorly  with  the  predictions  of  the  model 
as  shown  in  Fig.  5.12.  The  analytical  response  curve  follows  the  expected 
pattern  of  a  large  constant  stiffness  up  to  cracking,  followed  by  a  softer 
constant  stiffness  composed  of  the  addition  of  the  positive  stiffness 
contributed  by  the  steel  elements,  and  the  negative  stiffness  of  the  con¬ 
crete.  The  specimen  results  lack  the  initial  large  stiffness,  the  high 
cracking  load,  and  the  second  slope  of  smaller  stiffness.  The  greatest 
strength  contribution  by  the  concrete  is  only  approximately  40  kips, 
equivalent  to  80  psi,  whereas  the  model  predicts  the  concrete  will  carry 
up  to  100  kips,  and  the  concrete  strength  of  the  test  segment  is  exhausted 
at  a  small  strain,  leaving  only  the  steel  to  resist  the  load. 

5.1.3  General  Results 

The  results  for  the  specimens  with  prestressing  were  generally 
good.  The  model  response  curves  displayed  no  trend  in  being  either 
above  or  below  the  response  curves  of  the  test  specimens  and  therefore 
it  is  assumed  that  a  more  accurate  determination  of  the  effective  concrete 
tensile  strength  as  it  would  vary  between  specimens  would  provide  yet 


. 
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closer  agreement.  The  stiffness  of  both  the  model  and  the  test  specimen 
after  cracking  were  close,  indicating  that  the  postulated  declining 
branch  of  the  concrete  stress-strain  curve  in  tension  provides  results 
that  comply  with  the  test  observations. 

The  results  for  the  segment  with  an  applied  moment  were  comparable 
in  accuracy  with  the  results  for  the  segments  without  moments.  Based  on 
the  available  test  data,  the  predicted  moment-curvature  curves  would 
show  a  similarly  close  agreement. 

The  results  for  the  specimens  without  prestressing  were  poor,  al¬ 
though  the  sampling  may  be  small.  Compared  to  the  model  predictions, 
the  actual  specimens  were  soft.  Between  the  two  tests,  inconsistencies 
were  also  present.  One  exhibited  a  definite  cracking  load,  the  other 
did  not.  One  showed  the  concrete  to  persist  in  carrying  load  long  after 
"cracking",  while  the  other  displayed  a  rapid  decline  in  strength.  The 
presence  of  shrinkage  cracks,  which  would  not  be  a  factor  in  the  pre¬ 
stressed  specimens,  may  have  affected  the  strength  of  these  specimens. 
Nevertheless,  on  the  basis  of  these  two  tests,  it  appears  that  the 
effective  concrete  tensile  strength  for  reinforced  concrete  is  substan¬ 
tially  lower  than  that  for  prestressed  concrete. 

5.2  Comparison  of  Concrete  Constitutive  Relations 

The  test  wall  segments  were  also  analysed  using  the  uniaxial, 
bilinearly  elastic,  tension  cut-off  constitutive  relations  described  in 
Chapter  3.  The  results  of  these  analyses  are  also  shown  in  Figs.  5.4  - 
5.12.  In  addition,  the  diagrams  contain  the  contribution  of  the  steel, 
without  concrete,  including  prestressing,  towards  the  capacity  of  the 
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wall  segment  as  the  "steel  only"  curves. 

In  the  major  direction  of  Specimen  No.  1,  Fig.  5.4,  the  predicted 
cracking  load  from  the  tension  cut-off  analysis  is  the  same  as  that  from 
the  nonlinear  analysis,  and  both  are  above  the  observed  results.  After 
cracking,  the  concrete  no  longer  contributes  any  strength  by  the  tension 
cut-off  theory,  and  the  segment  response  by  this  theory  follows  the 
"steel-only"  solution,  which  is  lower  than  both  the  test  results  and  the 
nonlinear  results.  In  the  minor  direction,  the  predicted  cracking  load 
is  much  higher  than  the  results  from  both  the  nonlinear  analysis  and  the 
tests.  Again,  after  cracking,  the  tension  cut-off  solution  follows  the 
"steel-only"  curve,  which  is  below  the  test  and  nonlinear  results. 

The  predictions  from  the  tension  cut-off  analysis  are  high  in  one 
direction  of  Specimen  No.  3  (Fig.  5.5)  and  low  in  the  other.  In  the 
stiffer  direction,  the  cracking  load  is  lower  and  occurs  at  a  lower 
strain  than  the  test  and  nonlinear  results  show.  In  the  second  direction, 
the  cracking  load  is  higher  than  that  from  the  test  or  nonlinear  analysis. 

Specimen  No.  8  (Fig.  5.6)  is  similar  in  both  section  details  and 
loading  to  Specimen  No.  1.  The  comparison  between  the  bilinear,  nonlinear 
and  test  results  are  also  similar.  The  cracking  loads  of  the  two  analyses 

and  the  test  in  the  major  direction  are  close.  Again,  in  the  minor 

direction,  the  cracking  load  of  the  tension  cut-off  analysis  is  much 
higher  than  that  of  the  nonlinear  analysis  or  that  of  the  test  itself. 

The  predicted  cracking  loads  of  the  nonlinear  and  tension  cut-off 
analyses  are  the  same  for  Specimens  No.  5  and  6  shown  in  Figs.  5.7  and 

5.8.  For  Specimen  No.  5,  this  cracking  load  is  below  that  of  the  test 

and  for  Specimen  No.  6  the  predicted  cracking  load  is  above  that  of  the 


96 


test . 

For  the  segment  with  moment,  the  tension  cut-off  analysis  predicts 
the  correct  cracking  load  for  direction  1  on  Face  B  (Fig.  5.10)  and 
for  direction  2  on  Face  A  (Fig.  5.9).  After  cracking  the  results  of 
this  analysis  are  below  the  observed  results  for  both  curves,  as  most 
of  the  concrete  does  not  contribute  to  the  stiffness  of  the  segment 
after  this  point,  and  only  the  steel  remains  to  resist  the  loads.  Based 
on  Segment  Nos.  1,  3,  5,  6  and  8,  these  results  followed  the  expected 
pattern.  However  in  direction  1,  Face  A,  the  results  of  the  tension  cut¬ 
off  analysis,  after  cracking,  are  close  to  the  test  results,  implying 
that  the  concrete  of  the  test  segment  does  not  contribute  any  strength. 
Contradicting  this,  it  is  shown  from  Face  ^  direction  1  that  the  concrete 
is  stressed,  leading  to  the  conclusion  that  other  factors  may  have  been 
active  during  the  test  which  influenced  the  results.  In  direction  2  on 
Face  B,  the  tension  cut-off  results  closely  model  the  test  results. 

The  cracking  loads  for  Specimen  Nos.  4  and  7  (Figs.  5.11  and  5.12) 
are  again  the  same  from  both  analyses.  These  loads  are  significantly 
above  the  cracking  loads  observed  from  the  tests  on  these  specimens. 

By  the  bilinear  elastic,  tension  cut-off  constitutive  relations 
once  the  concrete  has  cracked  in  tension,  the  concrete  loses  all  its 
strength.  As  a  result,  in  the  tension  cut-off  analysis  of  each  specimen, 
once  the  segment  has  reached  its  cracking  load,  the  capacity  of  the  segment 
drops  off  to  the  "steel-only"  strength,  which  can  be  computed  in  a  very 
simple  manner.  The  predicted  amount  of  energy  which  the  segment  can 
absorb  is  therefore  lower  by  the  tension  cut-off  analysis  than  by  the 


nonlinear  analysis. 
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The  bilinear  elastic  constitutive  relations  formulated  here  are 
independent,  with  no  interaction  on  a  failure  curve.  Both  directions 
are  allowed  to  reach  the  full  tensile  cracking  load,  which  is  not 
usually  the  case  in  the  nonlinear  analysis.  If  the  segment  is  loaded 
in  only  one  direction,  then  the  cracking  load  of  both  the  nonlinear  and 
linear  analyses  are  comparable.  For  the  segments  loaded  in  both  direc¬ 
tions,  the  predicted  cracking  loads  in  the  direction  of  maximum  stress 
will  be  the  same  from  both  analyses.  These  comparable  results  will  be 
from  the  direction  where  the  concrete  stresses  at  failure,  in  the  non¬ 
linear  analysis  reach  the  full  tensile  strength.  In  the  orthogonal 
direction,  the  cracking  load  of  the  tension  cut-off  analysis  will  be 
greater  than  that  of  the  nonlinear  analysis.  So  in  general,  the  linear 
elastic,  uniaxial  tension  cut-off  analysis  will  predict  cracking  loads 
equal  to  or  greater  than  those  of  the  nonlinear  analysis.  However, 
exceptions  to  this  rule  can  occur  as  established  by  Specimen  No.  3. 
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TABLE  5.1  Summary  of  Wall  Segment  Data  -  Steel  Details 
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Specimen  Thickness  f'  f'  0.6  f'  Load 
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TABLE  5.2  Summary  of  Wall  Segment  Data  -  Concrete  Details 
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Ni 

<1>1 

n2 

m2 

(kips) 

(1/ in) 

(kips) 

(kip-in) 

84.0 

0.0 

42.0 

54.1 

382.0 

0.0 

141.0 

16  3.4 

335.0 

0.0 

267.5 

236.8 

535.0 

0.0 

300.0 

268.2 

535.0 

0.0 

340.0 

268.2 

TABLE  5.3 

Loading 

Sequence  for  Wall 

Segment  No.  12 
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Fig.  5.1  Typical  Test  Segment  Cross-Section  Orthogonal  to  Direction 


3)  7-6  Wire  Tendons 
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Fig.  5.2  Typical  Test  Segment  Cross-Section  Orthogonal  to  Direction 


Stress  in  Direction 
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Fig.  5.3 


Concrete  Stress  Path  for  Test  Segment  No.  1 
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Force-Deflection  Curve  of  Test  Segment  No. 
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Fig.  5.5  Force-Deflection  Curve  of  Test  Segment  No. 
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Fig.  5.6  Force-Deflection  Curve  of  Test  Segment  No. 
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Fig.  5.7  Force-Deflection  Curve  of  Test  Segment  No. 
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Fig.  5.8  Force-Deflection  Curve  of  Test  Segment  No. 
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Force-Deflection  Curve  of  Test  Segment  No.  12,  Face 
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Fig.  5.10  Force-Deflection  Curve  of  Test  Segment  No.  12,  Face 
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Fig.  5.11  Force-Deflection  Curve  of  Test  Segment  No. 
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Fig.  5.12  Force-Deflection  Curve  of  Test  Segment  No. 


6.  FURTHER  APPLICATIONS 


6.1  Introduction 

A  rigorous  and  widely  acceptable  method  of  analysing  a  concrete 
containment  vessel  would  involve  a  nonlinear  finite  element  or  finite- 
difference  technique.  However,  this  type  of  analysis,  if  used  as  a 
design  tool  can  be  costly  and  time  consuming.  In  this  chapter,  two 
simplified  alternative  design  checking  techniques  will  be  investigated, 
the  first  being  the  theory  and  resultant  program  developed  in  the  pre¬ 
vious  chapters  of  this  report,  which  considers  one  segment  of  the  total 
structure,  and  the  second  being  the  interaction  diagrams  presented  by 
Epstein  and  Murray  (1976) ,  which  considers  only  one  direction  of  one 
segment.  These  two  methods  will  then  be  compared  with  the  results 
obtained  from  a  B0S0R5  analysis,  which  considers  the  entire  structure 
and  for  which  the  program  was  modified  to  include  an  elastic-plastic 
constitutive  model  for  concrete.  Before  making  these  comparisons,  an 
explanation  of  the  construction  and  interpretation  of  the  interaction 
plots  will  be  given. 

6.2  Interpretation  of  Interaction  Diagrams 

As  an  alternative  to  the  type  of  analysis  discussed  in  this  report, 
in  which  a  segment  is  subjected  to  a  prescribed  loading  sequence,  Epstein 
and  Murray  (1976)  constructed  interaction  diagrams.  From  those  diagrams, 
plotted  in  axial  force-moment  space,  the  various  conditions  of  first  crack¬ 
ing,  first  yield,  through  cracking,  second  yield,  tendon  rupture  and 
crushing  of  the  concrete  can  be  predicted,  not  just  for  a  single  combina¬ 
tion  of  loads,  but  for  any  combination.  Figure  6.1  contains  a  typical 
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interaction  diagram  from  that  report. 

_ In  constructing  the  interaction  plots,  a  nonlinear  one-dimensional 
concrete  constitutive  relation  was  employed.  Each  plot  ignores  biaxial 
effects  and  is  restricted  to  a  given  set  of  parameters,  namely:  the 
concrete  strength,  steel  strength  and  the  segment  geometry.  In  addition, 
no  thermal  or  creep  conditions  could  be  considered,  and  it  was  assumed 
that  no  prior  cracking  had  occurred  due  to  thermal  or  shrinkage  stresses. 

To  use 'the  interaction  plots  the  effects  of  prestressing  must  be 
included  as  an  external  load,  even  though  the  strains  in  the  tendons 
resulting  from  prestressing  are  implicitly  considered  in  each  diagram. 
From  the  location  representing  the  axial  load  and  moment  induced  by  the 
prestressing,  the  live  loads  are  added.  Due  to  the  nature  of  the  plot, 
the  live  loads  must  be  in  the  form  of  axial  loads  and  moments,  not  in 
the  form  of  imposed  strains.  Predictions  from  the  interaction  diagrams 
can  be  made  for  a  load  path  in  which  the  ratio  of  axial  load  to  moment  is 
kept  constant.  Since  the  actual  loading  history  in  a  structure  is  a 
curved  load  path,  especially  after  first  cracking,  the  straight  line 
assumption  may  lead  to  erroneous  results. 

As  presented  by  Epstein  and  Murray  (1976),  the  interaction  diagrams 
may  create  some  confusion.  This  can  be  illustrated  by  constructing  the 
load  paths  Nos.  1  and  2  as  shown  in  Fig.  6.1.  If  the  segment  loading 
follows  path  No.  1,  then  the  predictions  of  the  plots  are  clear.  The 
prestressing  brings  the  loading  condition  to  point  A.  As  the  live  load 
is  applied,  the  segment  begins  to  crack  at  B.  Loading  continues  and  at 
point  C,  one  of  the  layers  of  reinforcing  steel  begins  to  yield.  Failure 
of  the  segment  occurs  at  point  D  when  crushing  of  the  concrete  is  first 
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experienced  on  the  outside  layer. 

-When  following  load  path  No.  2  however,  the  nature  of  the  response 
is  not  so  evident.  Again  beginning  at  point  A,  the  first  limit  state 
encountered  is  through  cracking  at  point  E,  indicating  that  open  cracks 
exist  completely  through  the  section.  As  the  loading  continues  from 
point  E,  first  cracking  is  observed  at  point  F.  First  yield  occurs 
at  point  G,  through  cracking  is  again  indicated  at  point  H,  second  yield 
at  point  I  and  final  failure  results  from  tendon  rupture  at  point  J. 

This  apparent  sequence,  in  which  first  cracking  occurs  after  open 
through  cracking,  is  obviously  incorrect,  and  a  proper  interpretation 
requires  a  consideration  of  the  load  history. 

The  area  of  the  diagram  bounded  by  the  through  cracking  envelop  is 
the  only  location  where  load  history  must  be  considered.  If  the  section 
had  not  yet  been  subject  to  any  loads,  the  following  events  would  occur 
along  load  path  No.  2.  The  first  limit  state  to  be  encountered  after 
leaving  A  would  again  be  open  through  cracking,  but  now  at  location  F. 

As  the  loading  continues,  first  yield  would  be  observed  at  point  G,  and 
open  through  cracking  again  at  point  H.  At  point  H,  the  segment  no 
longer  contains  open  through  cracks  as  the  moment  has  forced  one  side 
of  the  segment  into  compression.  This  condition  however,  is  accompanied 
by  large  strain  gradients.  The  lines  labelled  through  cracking  form 
an  envelop  within  which  the  condition  of  open  through  cracks  exist. 
Therefore  at  location  H,  the  loading  leaves  the  open  through  cracking 
envelop  and  the  cracks  on  one  side  of  the  segment  close.  After  point  H, 
second  yield  occurs  at  point  I  and  as  before,  failure  is  the  result  of 
tendon  rupture  at  point  J. 
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Load  reversal,  however,  will  lead  to  different  results.  Once  the 
loading  has  brought  the  forces  and  moments  to  some  location  after  point 
F,  where  first  cracking  and  open  through  cracking  occur  simultaneously, 
if  some  of  the  load  is  removed  the  cracks  will  not  close  until  the  through 
cracking  envelop  is  reached  at  point  E.  If  the  segment  has  not  previously 
been  cracked,  then  first  cracking  and  open  through  cracking  occur 
simultaneously  at  F.  However,  once  the  segment  has  cracked,  the  cracks 
will  not  close  until  E  or  H  is  reached  and  with  reloading,  the  open 
cracks  will  appear  again  at  E. 

Within  the  open  through  cracking  envelop,  only  the  two  layers  of 
reinforcing  steel  and  the  prestressing  steel  are  resisting  the  axial 
forces  and  moments.  This  results  in  three  other  characteristics  of 
the  interaction  diagrams.  First,  the  vertical  line  which  extends  from 
the  lower  corner  of  the  through  cracking  envelop  to  the  tendon  rupture 
line  is  part  of  the  through  cracking  envelop.  If  the  load  combination  is 
on  this  vertical  line,  then  open  through  cracks  exist,  and  both  layers 
of  steel  have  yielded.  Secondly,  the  change  in  slope  of  the  through 
cracking  line  below  the  first  yield  line  results  from  yielding  of  the 
prestressing  strands.  The  stress-strain  law  considered  for  the  pre¬ 
stressing  is  linear  up  to  0 . 7  a  .  After  this,  the  tangential  modulus 

pu 

declines  to  zero  at  a  .  Therefore,  at  the  "kink"  in  the  through  cracking 

pu 

line,  the  prestressing  strands  have  reached  the  yield  stress.  Third,  both 
of  the  steel  layers  have  yielded  on  the  vertical  line  extending  from  the 
junctions  of  the  first  yield  lines  to  the  second  yield  lines,  and  this 
line  should  also  be  labelled  second  yield.  As  this  example  section  has 


steel  layers  of  unequal  area,  this  line  is  not  on  the  moment  axis  (M=0) 
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The  interaction  diagrams  presented  by  Epstein  and  Murray  (1976) 
are  correct  but  must  be  used  with  care,  especially  in  the  region  bounded 
by  the  through  cracking  lines.  If  the  section  is  previously  uncracked, 
cracking  will  not  occur  anywhere  on  the  diagram  until  the  first  cracking 
line  is  reached.  If  this  line  is  reached  within  the  through  cracking 
envelop,  then  first  cracking  and  through  cracking  occur  simultaneously. 

If  the  segment  has  already  been  cracked,  then  the  concrete  has  no  strength 
in  tension  and  through  cracks  will  exist  if  the  load  state  is  anywhere 
within  the  through  cracking  envelop. 

6.3  Analysis  of  Containment  Structure  Segments 

To  further  assess  the  abilities  of  the  segment  analysis  as  developed 
herein,  comparisons  were  made  indirectly  with  a  scale  model  test  structure 
which  was  tested  at  the  University  of  Alberta  in  1978.  This  test  struc¬ 
ture  was  first  analysed  using  30S0R5  which  provided  an  energy  finite 
difference  solution  of  the  axisymmetric  containment  vessel  subject  to 
internal  pressure.  For  this  analysis  of  the  test  structure,  one  slice 
of  the  structure,  composed  of  102  segments,  was  examined  (Simmonds,  et  al, 
1979).  The  results  of  the  B0S0R5  analysis  compare  well  with  the  test 
structure  and  shall  therefore  be  assumed  to  correctly  represent  the  test 
model.  The  segment  analysis  technique  of  this  report  is  now  applied  to 
segments  from  the  test  structure  and  the  results  compared  to  the  B0S0R5 
analysis . 

For  the  comparison,  four  segments  were  taken  from  the  B0S0R5 
analysis.  These  segments  are  labelled  9-10,  3-2,  3-16  and  7-4.  (The 
notation,  of  the  form  m-n,  arises  from  the  B0S0R5  model,  in  which  m  is 


- 


118 


the  shell  component,  and  n  is  the  station  within  the  component.) 

Segment  9-10  is  at  the  top  of  the  dome,  Segment  3-2  is  at  mid  height  in 
the  cylinder  wall,  Segment  3-16  is  at  the  junction  of  the  cylinder  with 
the  ring-beam  and  Section  7-4  forms  the  junction  of  the  ring-beam  with 
the  dome.  The  details  of  the  segment  geometries  are  listed  in  Tables 

6.1  and  6.2,  along  with  pertinent  material  parameters. 

The  results  of  the  segment  analyses  are  presented  such  that  the 
forces  and  moments  on  orthogonal  load-carrying  faces  are  examined  in¬ 
dividually.  These  faces  are  designated  as  meridional  (or  vertical)  by 
M,  and  circumferential  by  C.  That  is,  the  forces  on  face  C  act  in  the 
circumferential  direction  and  face  C  is  orthogonal  to  these  forces. 

For  Segment  9-10,  at  the  top  of  the  dome,  the  segment  geometry  and  loads 
are  similar  for  both  faces  and  therefore  only  the  meridional  face  is 
discussed.  The  results  of  the  circumferential  face  are  similar. 

The  graphical  results  of  the  various  solutions  are  shown  in  Figs. 

6.2  -  6.8.  These  figures  deviate  from  the  sign  convention  adhered  to 

in  previous  pages  herein  in  that  compressive  axial  forces  are  considered 
positive.  The  solid  lines  which  form  the  interaction  diagrams  represent 
the  uniaxial  elastic  solutions  of  the  Epstein  and  Murray  report  (1976). 
The  solutions  of  the  two-dimensional  nonlinear  analysis,  and  of  the 
B0S0R5  nonlinear  analysis  are  shown  as  discrete  points  plotted  on  the 
interaction  diagram.  Each  point  represents  an  increment  of  10  psig  of 
internal  pressure  within  the  vessel  starting  from  the  reference  state 
which  includes  prestressing,  gravity,  and  the  weight  of  the  water  which 
fills  the  structure.  The  points  for  the  B0S0R5  analysis  are  plotted  up 
to  130  psig.  By  the  B0S0R5  analysis,  the  predicted  failure  load  occurs 
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at  approximately  134  psig. 

-While  the  live  loading  required  for  the  B0S0R5  solution  is  specified 
as  an  internal  pressure,  the  wall  segment  analysis  developed  in  this 
report  requires  a  loading  history  in  terms  of  axial  forces  or  strains, 
moments  or  curvatures.  The  source  of  this  loading  path  must  come  from 
a  solution  which  considers  the  entire  structure,  not  just  its  elements  in 
isolation.  For  the  following  comparisons,  the  results  of  the  B0S0R5 
analysis  up  to  an  internal  pressure  of  30  psig,  at  which  point  all  of 
the  materials  were  behaving  elastically,  were  used  to  define  the  loading 
history  required  by  the  wall  segment  analysis.  Two  different  loading 
histories  were  used.  The  first  considers  a  constant  ratio  of  the  moments 
to  axial  forces  on  a  face  which  result  from  the  internal  pressures,  and 
these  loads  are  added  to  the  reference  state.  This  shall  be  called  the 
N-M  analysis.  The  second  loading  history  considers  that  the  axial  forces 
and  curvatures,  resulting  from  the  internal  pressure  have  a  constant 
ratio.  This  method  shall  be  called  the  N-cf>  analysis.  For  both  of  these 
loading  histories,  the  ratio  of  the  change  in  axial  force  to  the  change 
in  internal  pressure  was  assumed  to  be  constant  for  the  entire  loading 
sequence . 

The  load  path  for  a  direct  analysis  by  the  interaction  plots  must 
also  be  specified  in  terms  of  axial  forces  and  moments.  Here  the  axial 
force  to  moment  ratio  resulting  from  the  elastic  response  will  be  extra¬ 
polated  until  a  failure  condition  is  reached  on  the  interaction  plot, 
and  knowing  the  increase  in  loads  for  a  given  increase  in  internal 
pressure,  the  pressure  at  which  this  failure  condition  has  been  reached 
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will  be  estimated.  The  results  of  this  and  the  other  procedures  are 
listed  in  Table  6.3.  This  table  lists  the  maximum  possible  internal 
pressure  that  can  be  reached  by  each  analysis,  along  with  the  axial 
force  and  moment  at  this  level  of  internal  pressure.  Since  the  B0S0R5 
analysis  deals  with  the  entire  structure,  the  maximum  load  to  which  any 
segment  can  be  subjected  is  the  load  at  which  the  weakest  segment  fails. 
The  B0S0R5  analysis  failed  to  converge  at  pressures  above  134  psi. 

B0S0R5  values  shown  on  Table  6.3  are  extrapolated  pressures  required  to 
produce  failure  of  the  tendons  at  the  segment  location.  For  the  N-M 
and  N— 4>  analyses,  the  maximum  load  to  which  any  segment  can  be  subjected 
is  that  load  at  which  the  weaker  direction  of  the  segment  fails.  Since 
the  interaction  plots  are  uniaxial,  each  direction  of  a  segment  fails 
at  different  load  levels.  If  the  smallest  ultimate  load  of  the  four 
segments  is  taken  from  each  of  the  analyses,  the  maximum  loads  as  pre¬ 
dicted  by  the  B0S0R5,  N-<j>,  N-M  and  interaction  analyses  are  138,  120,  95 
and  100  psig. 

6.3.1  Segment  9-10 

The  results  from  segment  9-10  are  shown  in  Fig.  6.2.  Since  the 
segment  details  and  loads  are  the  same  in  both  directions,  this  figure 
represents  the  results  for  both  load-carrying  faces  of  segment  9-10. 

The  solid  lines  represent  the  familiar  interaction  diagram  while  the 
other  analyses  are  shown  by  their  designated  symbols.  All  the  segments 
have  been  analysed  assuming  the  dimensions  are  one  inch  wide  by  one  inch 
long,  and  this  is  reflected  in  the  units  of  the  axes.  The  thicknesses 
of  the  segments  vary  and  are  listed  in  Table  6.1. 
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The  straight  line  load  paths  of  the  N-M  segment  analyses  and  the 
constant  ratio  extrapolations  on  the  interaction  plots  will  always  have 
the  same  direction,  and  in  this  instance  both  indicate  failure  at  an 
internal  pressure  of  110  psi.  This  agreement  in  results  will  not  always 
occur  but  it  is  more  likely  to  happen  when  the  orthogonal  load-carrying 
faces  of  the  segment  have  similar  loadings  and  details.  However,  the 
solutions  obtained  by  assuming  a  constant  ratio  of  axial  force  to  moment 
do  not  follow  the  B0S0R5  solution  as  well  as  the  N-<J>  analysis.  After 
crossing  the  first  cracking  line  of  the  interaction  diagram,  both  the 
N— 4>  and  B0S0R5  results  predict  that  the  moments  decrease  to  almost 
zero.  However,  the  N— <f>  analysis,  which  assumes  that  the  axial  force 
increases  linearly  with  the  internal  pressure  predicts  failure  of  the 
segment  at  120  psig  whereas  the  B0S0R5  solution,  which  is  not  restricted 
by  this  assumption,  indicates  that  the  segment  can  withstand  at  least 
138  psig.  The  assumption  of  linearity  between  membrane  force  and 
pressure  is  reasonably  accurate  in  this  instance  but  it  is  not  sufficiently 
precise  to  make  the  agreement  in  internal  pressures  at  failure  between 
these  two  methods  any  closer. 

6.3.2  Segment  3-2 

The  results  of  the  forces  and  moments  in  the  meridional  direction  of 
segment  3-2  are  shown  in  Fig.  6.3  and  the  circumferential  direction  is 
shown  in  Fig.  6.4.  Segment  3-2  is  located  half-way  up  the  cylindrical  wall 
of  the  structure  and  is  little  affected  by  the  base  or  the  ring-beam  at 
the  top.  In  the  meridional  direction  the  section  is  only  subjected  to 
tensile  forces;  the  moments  are  very  small.  As  a  result  all  of  the 


. 
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analyses  have  similar  predictions.  The  points  of  the  B0S0R5,  N-M  and  N— <f> 
analysis  do  not  reach  the  outer  edge  of  the  interaction  plot  as  the  failure 
occurs  first  in  the  orthogonal  circumferential  direction.  In  this  direc¬ 
tion,  the  N-M  analysis  again  has  moments  close  to  zero.  However,  after 
cracking  both  the  N— <f>  and  B0S0R5  analyses  predict  increasing  moments 
across  this  section  which  is  eccentrically  loaded  by  the  prestressing. 
Failure  occurs  in  this  direction  when  the  prestressing  tendons  fracture. 

For  this  segment,  as  in  segment  9-10,  the  closest  agreement  with 
the  B0S0R5  results  is  obtained  from  the  N-cj>  analysis.  Also,  in  common 
with  segment  9-10,  this  segment  is  not  located  close  to  any  disconti¬ 
nuities  . 

6.3.3  Segment  3-16 

Segment  3-16  is  located  at  the  top  of  the  cylindrical  wall 
immediately  next  to  the  top  ring  beam.  The  forces  and  moments  applied 
to  the  segment  in  the  meridional  direction  are  shown  in  Fig.  6.5  and  in 
the  circumferential  direction  in  Fig.  6.6.  Arising  from  the  proximity 
of  this  segment  to  the  ring  beam,  large  moments  must  be  resisted  in  the 
meridional  direction.  Whereas  the  N—  analysis  indicates  that  the 
moments  are  reduced  after  cracking,  the  B0S0R5  analysis  shows  that  these 
moments  must  be  carried  past  the  initiation  of  cracking  and  beyond  the 
point  of  first  yielding.  For  this  section,  the  N-M  analysis  provides 
better  agreement  with  the  B0S0R5  results  than  does  the  N— 4>  analysis. 

In  the  circumferential  direction,  the  comparatively  large  and  con¬ 
sequently  stiff  ring  beam  absorbs  the  majority  of  the  forces  and  therefore 
the  results  are  centred  near  the  origin  of  the  interaction  plot,  away 
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from  any  failure  modes.  Here  again  though,  like  the  meridional  direction, 
the  N-M  analysis  provides  results  which  are  closer  to  the  B0S0R5  results 
than  the  results  of  the  N-cj>  analysis. 

6.3.4  Segment  7-4 

This  segment  is  located  near  the  bottom  of  the  dome  which  is  linked 
to  the  cylinder  through  the  ring  beam.  This  segment  therefore  also 
reflects  many  of  the  effects  that  were  seen  in  segment  3-16.  The  results 
for  the  meridional  direction,  shown  in  Fig.  6.7  parallel  those  in  the 
meridional  direction  of  segment  3-16.  After  first  cracking,  the  moments 
predicted  by  the  N-<J>  analysis  drop  off.  It  is  apparent  though,  that 
these  moments  must  be  resisted,  as  the  B0S0R5  analysis  shows  their  presence 
after  the  initiation  of  cracking,  after  first  yielding,  up  to  crushing. 

This  behavior  is  closer  to  the  response  demonstrated  by  the  N-M  assumptions 
than  by  the  N-<f>  assumptions . 

The  results  in  the  circumferential  direction,  shown  in  Fig.  6.8  are 
also  similar  to  those  of  section  3-16C.  The  large  adjoining  ring  beam 
attracts  a  large  proportion  of  the  forces  during  the  initial  stages  of 
loading.  However,  the  B0S0R5  results  do  not  adhere  to  the  assumption 
made  for  the  N-M  and  N-<f>  analyses  that  the  axial  force  is  proportional 
to  the  internal  pressure.  As  shown  by  the  B0S0R5  points  plotted  in 
Fig.  6.8,  which  represent  increments  of  10  psig  of  internal  pressure, 
the  change  in  axial  force  for  a  constant  change  in  internal  pressure 
increases  at  higher  pressures.  As  a  result  of  this  assumption,  changing 
ratios  of  axial  force  to  internal  pressure  cannot  be  predicted. 
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6.3.5  Summary  of  Results 

-Two  different  types  of  behavior  were  observed.  When  the  segment 
investigated  is  removed  from  any  discontinuities,  then  the  N-cf>  analysis 
gave  results  which  most  closely  approximate  the  B0S0R5  results.  However, 
when  the  segment  was  close  to  a  discontinuity,  then  the  N-M  analysis  was 
the  more  accurate.  For  the  most  part,  the  B0S0R5  results  were  bounded  by 
the  N-M  and  N~4>  analyses  on  the  interaction  curve.  The  N-M  analysis  gave 
an  underestimate  for  ultimate  load  in  every  instance  whereas  the  N-<J> 
analysis  could  predict  higher  or  lower  ultimate  loads.  The  results 
obtainable  from  the  uniaxial  nonlinear  elastic  interaction  charts  are 
nearly  identical  to  those  of  the  biaxial  nonlinear  N-M  analyses. 
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Fig.  6.1  Typical  Interaction  Diagram 
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Fig.  6.2  Test  Structure  Segment  9-10 
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Fig.  6.4  Test  Structure  Segment  3-2C 
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Fig.  6.5  Test  Structure  Segment  3-16M 
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Fig.  6.6  Test  Structure  Segment  3-16C 
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Fig.  6.7  Test  Structure  Segment  7-4M 


135 


(•m/sdpi)  aojoj[  jeixy 


Fig.  6.8  Test  Structure  Segment  7-4C 


7.  SUMMARY  AND  CONCLUSIONS 


7.1  "Summary 

A  computer  model  was  formulated  to  analyze  a  concrete  wall  segment 
containing  two  degrees  of  freedom  in  the  direction  of  each  of  the  two 
principal  curvatures.  In  addition  to  the  concrete  matrix,  the  segment 
may  include  uniaxial  elements  of  reinforcing  bars  and  prestressing 
strands . 

The  segment  analysis  can  be  accomplished  with  either  one  of  two 
available  concrete  constitutive  relations.  The  first  is  a  bilinear 
elastic  uniaxial  relation  incorporating  a  two-part  stress-strain  curve 
in  the  compression  region,  and  a  tension  cut-off.  The  second  concrete 
constitutive  relation  is  a  biaxial  nonlinear  elastic  type.  In  this 
relation,  the  real  strains  are  first  converted  to  equivalent  uniaxial 
strains  (Darwin  and  Pecknold,  1974).  Then  for  a  given  stress  ratio  the 
corresponding  point  on  the  failure  curve  is  found  (Kupfer  and  Gerstle, 
1973),  the  appropriate  value  for  Poisson’s  ratio  is  determined  (Elwi 
and  Murray,  1979),  and  from  curves  developed  in  this  report,  the  equiva¬ 
lent  uniaxial  strains  at  ultimate  stress  are  ascertained.  The  current 
stresses  are  then  computed  from  a  stress-strain  equation  which  includes 
a  softening  branch  (Saenz,  1964) . 

The  constitutive  relation  of  the  reinforcing  steel  consists  of  the 
conventional  two-part  elas tic-perf ectly  plastic  curve  whereas  the  res¬ 
ponse  of  the  prestressing  steel  is  described  by  a  piece-wise  linear 
stress-strain  curve. 

The  segment  can  be  subjected  to  four  different  types  of  load. 
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These  loads  include  prestressing  forces,  creep  effects,  thermal  effects 
and  live  loads.  The  solutions  under  these  loads  are  found  by  the  Newton- 
Raphson  technique.  The  response  of  the  segment  to  a  given  set  of  strains 
is  found  by  layering  the  segment,  giving  each  layer  appropriate  properties 
and  then  integrating  the  responses  from  the  individual  layers  to  find  the 
segment  response.  Corrections  are  made  to  this  set  of  strains  and  the 
process  repeated. 

Comparisons  were  then  made  between  this  mathematical  model  and 
results  from  wall  segment  tests.  Both  the  bilinear  and  nonlinear  solu¬ 
tions  were  compared  with  seven  test  segments  involving  combinations  of 
pure  axial  loads,  and  one  test  segment  subject  to  axial  loads  and  a 
moment . 

The  model  was  then  applied  to  predict  the  nonlinear  response  of 
wall  segments  from  a  containment  vessel  test  structure.  Using  the 
results  of  the  B0S0R5  computer  program,  the  initial  axial  force  to  moment 
ratios  and  the  initial  axial  force  to  curvature  ratios  were  determined. 

It  was  then  assumed  that  either  of  these  relations  would  remain  constant 
during  the  entire  loading  sequence.  These  assumptions  formed  the  basis 
of  the  N-M  and  N-<J>  analyses  of  Chapter  6.  The  results  from  these 
simplified  analyses  were  then  compared  with  those  from  a  nonlinear 
B0S0R5  analysis  of  the  entire  structure. 

7.2  Conclusions 

The  mathematical  model,  complete  with  the  nonlinear  biaxial  con¬ 
crete  constitutive  relation,  which  has  been  presented  in  this  report  to 
analyse  reinforced  and  prestressed  wall  segments  will  predict  with 
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acceptable  accuracy  the  response  of  a  prestressed  wall  segment  to 
biaxial  membrane  forces,  biaxial  moments,  or  both.  The  predicted  res¬ 
ponse  of  reinforced  wall  segments,  as  opposed  to  prestressed  wall  segments 
is  not  close  to  the  test  results.  However  these  discrepancies  may  be  the 
result  of  a  conflict  between  the  method  of  loading  and  the  assumptions 
of  the  segment  loading  analysis.  More  importantly,  it  is  probable 
that  there  are  shrinkage  cracks  and  stresses  in  the  reinforced  segments 
which  were  not  present  in  the  prestressed  segments. 

The  wall  segment  model  used  with  the  bilinear  elastic,  tension  cut¬ 
off  concrete  constitutive  relation  provides  results  which  may  or  may  not 
be  accurate  up  to  the  cracking  load,  and  which  underestimate  the  post 
cracking  stiffness  of  the  wall  segment  in  all  cases. 

The  models  were  not  verified  in  the  compressive  failure  zones  and 
no  conclusions  can  be  drawn  as  to  their  performance  under  compressive 
loading . 

Although  the  wall  segment  analyses  labelled  N-M  and  N-tj)  bracket  the 
results  obtained  from  a  nonlinear  finite  element  program  on  an  interaction 
diagram,  neither  analysis  provides  good  agreement  for  all  segments  which 
might  be  extracted  from  the  structure.  For  the  wall  segments  unhindered 
by  boundary  elements,  the  N— analysis  generates  the  best  results.  For 
the  segments  closer  to  the  boundary  elements,  the  N-M  analysis  provides 
better  results.  In  all  cases,  the  N-M  analysis  provides  an  underestimate 


of  the  maximum  strength. 
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APPENDIX  A 


COWSAC  USER'S  MANUAL 


A.l  Program  Objectives 

This  program  is  designed  to  predict  the  biaxial  response  of  an 
arbitrarily  prestressed  and  steel  reinforced  concrete  wall  segment  of  any 
dimension,  in  a  plane  stress  condition,  including  out-of-plane  curvatures 
with  respect  to  creep,  temperature  and  live  loads.  Given  a  force  or 
moment  on  any  face  orthogonal  to  the  principal  radii  of  curvature,  the 
program  will  find  the  resulting  strains  and  curvatures,  or  vice  versa. 
Thus,  inputting  any  combination  of  forces  and/or  displacements  on  the 
two  faces  of  the  segment  will  result  in  the  solution  of  the  unknown 
displacements  and/or  forces. 

A. 2  Input  Deck 

Input  is  to  be  read  from  unit  5  and  shall  conform  to  the  description 
given  below.  With  the  MTS  system,  data  can  be  free-format  with  commas 
separating  entries,  or  shall  follow  the  specific  format  given.  All 
input  quantities  should  be  in  consistant  units. 


A. 2.1  Data  Deck 

Description 

NUMBER  OF  CARDS 
(or  lines) 

FORMAT 

VARIABLES 

1 

20A4 

TITLE 

1 

13,  5F10.4 

NLAYER,  T,  WIDTH,  HEIGHT, 
Rl,  R2 

1 

13,  2E12.5 

NREBAR,  ER,  FYR 

NREBAR 

13,  2F10.4 

NRD(I)  ,  AREAR(I)  ,  RCEN(I) 
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NUMBER  OF  CARDS 

FORMAT 

VARIABLES 

1 

15,  4E12.7 

NFLAGD ,  ALOADD(l)  ,  ALOADD(2)  , 
AL0ADD(3) ,  ALOADD(4) 

1 

13,  2E10.5, 

413,  F10.5 

NPRE,  EP,  FYP ,  KP,  NPSTEP, 

KDO,  NCRP ,  FCRP 

NPRE 

13,  3F10.4, 

13 

NPD(J) ,  AREAP(J) ,  PCEN(J) , 
PFORCE(J),  NPN(J) 

NPSTEP 

15,  4E12.7 

NFLAGL(I) ,  ALOADL(I,l) ,  ALOAD(I,2) 
ALOAD(I , 3) ,  ALOAD(I , 4) 

1 

514 

NLOADS ,  NJUMPX,  NPDO,  NTHERM, 
NEPSTE 

1,  if  NTHERM 
^  0 

5E12.5,  14 

Tl,  T2,  T3 ,  ALPHAC ,  ALPHAS, 

NTCODE 

1 

15,  4E12.5, 

15 

NFLAG(I) ,  ALOAD(I , 1) ,  ALOAD(I,2), 
ALOAD(I , 3) ,  ALOAD(I , 4) ,  NINC(I) 

the  next  two  cards  are  to  be  repeated  NLOAD  -  1  times 


1 

14 

NJUMPX 

1 

15,  4E12.5,  15 

NFLAG(I) ,  ALOAD(I , 1) ,  AL0AD(I,2), 
ALOAD(I , 3) ,  ALOAD(I , 4) ,  NINC(I) 

the 

last  card  shall  always  be 

1 

4F10.4,  215,  2E12.7 

FC,  FTC,  POISS,  CEMOD,  NTYPE , 

NSUB ,  EIO ,  EIOT 

A. 2. 

2  Explanation  of  Variables 

1. 

TITLE 

any  alphanumeric  string  of  80  or  less 

characters . 

2. 

NLAYER 

number  of  concrete  layers  into  which 
segment  thickness  is  to  be  subdivided 

the  wall 

3. 

T 

wall  segment  thickness. 

4. 

WIDTH 

dimension  of  wall  segment  parallel  to 

direction  2. 

5. 

HEIGHT 

dimension  of  wall  segment  parallel  to 

direction  1  . 

’ 
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6.  Rl 

radius  curvature  of  wall  about  a  line  parallel  to 
direction  1.  If  no  curvature  is  desired  (i.e. 
a  flat  wall)  give  Rl  a  negative  value. 

7.  R2 

radius  curvature  of  wall  about  a  line  parallel 
to  direction  2.  If  no  curvature  is  desired 
give  R2  a  negative  value. 

8 .  N REBAR 

total  number  of  reinforcing  steel  layers. 

9 .  ER 

elastic  modulus  of  reinforcing  steel.  If  NREBAR 
equals  0,  ER  and  EYR  may  be  omitted. 

10.  FYR 

yield  stress  of  reinforcing  steel. 

11.  NRD(I) 

either  1  or  2;  specifies  the  direction  parallel 
to  the  steel  bars,  for  layer  I. 

12.  AREAR(I) 

the  total  area  of  steel  bars  in  layer  I. 

13.  RCEN(I) 

the  distance  from  the  centre  line  of  the  concrete 
wall  to  the  middle  of  the  steel  layer,  a  positive 
distance  measured  moving  away  from  the  centre  of 
curvature . 

14 .  NFLAGD 

load  code  for  dead  load  (see  Section  A. 2. 2  entitled 
Load  Codes).  To  NFLAGD  only,  add  1000. 

15.  ALOADD(J) 

the  4  dead  load  forces  or  strains  to  be  put  on 
the  wall  segment.  The  loads  shall  be  forces  or 
strains  as  specified  by  NFLAGD. 

16.  NPRE 

the  number  of  prestressing  layers. 

17.  EP 

elastic  modulus  of  the  prestressing  strands. 

18.  FYP 

the  ultimate  strength  of  the  prestressing  strands. 
From  this  and  EP ,  a  stress-strain  curve  for  the 
strands  is  calculated. 

19 .  KP 

either  zero  (0)  or  one  (1).  If  KP  =0,  then  the 
input  prestressing  force  is  before  elastic  res¬ 
ponse.  If  KP  =  1,  then  the  input  prestressing 
force  is  after  elastic  response. 

20.  NPSTEP 

the  number  of  steps  required  to  apply  prestressing. 
If  five  tendons  are  released  individually,  not 
simultaneously,  then  there  would  be  five  pre¬ 
stressing  steps. 

21.  KDO 

if  KDO  =  999,  do  not  apply  prestress. 

145 


22.  NCRP 

if  NCRP  is  not  equal  to  zero,  then  apply  creep 
under  constant  load  after  prestressing  steps 
are  complete. 

23.  FCRP 

creep  factor  FCRP  -  E.  /  E  ,  _  _ 

long  term  short  term 

24.  NPD(J) 

either  1  or  2;  specifies  the  direction  which  the 
prestressing  strands  parallel,  for  layer  J. 

25.  A REAP (J) 

the  total  area  of  prestressing  strands  in  layer  J 

26.  PCEN(J) 

the  distance  from  the  centre  line  of  the  concrete 
wall  to  the  middle  of  the  strand  layer,  positive 
away  from  the  centre  of  curvature. 

27.  PFORCE(J) 

prestressing  force  applied  to  layer  J  consistent 
with  the  value  of  KP  above. 

28.  NPN(J) 

the  position  in  the  prestressing  sequence  in 
which  prestressing  is  applied  to  the  layer  J, 
i.e.  if  NPN(J)  is  4,  then  the  prestressing  force 
is  applied  to  this  layer  in  the  fourth  step 
out  of  NPSTEP  steps. 

29.  NFLAGL(I) 

load  code  for  prestressing  loads  AL0AD(I,J). 

30.  ALOADL ( I , J ) 

the  load  vector  to  which  the  segment  is  subject 
after  each  prestress  step  I. 

31.  NLOADS 

the  number  of  live  loads. 

32.  NJUMPX 

the  number  of  increments  into  which  the  following 
live  load  is  to  be  divided. 

33.  NPDO 

if  NPDO  is  set  equal  to  0,  then  a  summary  of  the 
stresses  in  each  layer  is  output  after  each  load 
increment . 

34 .  NTHERM 

if  set  equal  to  1,  2,  or  3,  then  thermal  loads 
are  applied  to  the  segment  before  the  live  loads. 
The  temperature  distribution  can  be  one  of  three: 

NTHERM  =1,  AT  =  T1  +  T2*  Z  +  T3-Z2 

T3 

NTHERM  =  2,  AT  =  T1  +  T2*Z 

NTHERM  =  3,  AT  =  T1  +  T2*eZ‘T3 

For  these  equations  only,  the  distance  coordinate 
Z  is  defined  as  zero  on  the  inside  surface, 
positive  moving  towards  the  outside. 
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35.  NEPSTE 

if  NEPSTE  =  1,  then  a  ratio  of  axial  force  to 
moment  should  be  input  where  curvatures  are 
expected  in  the  load  vector.  The  curvature  for 
a  unit  applied  moment  is  calculated,  and  the 
force-moment  ratio  is  converted  to  a  force- 
curvature  ratio.  From  this  ratio  an  applied 
curvature  for  a  given  force  is  calculated  and 
the  load  vector  of  axial  forces  and  curvatures 
is  constructed. 

36.  Tl,  T2 ,  T3 

thermal  load  constants  applicable  in  the  equations 

to 

above.  The  entire  card  containing  Tl  through 
NTCDDE  must  be  omitted  if  thermal  loads  are  not 

38. 

desired . 

39.  ALPHAC 

thermal  coefficient  of  expansion  for  concrete. 

40.  ALPHAS 

thermal  coefficient  of  expansion  for  steel. 

41.  NTCODE 

load  code  for  thermal  loads. 

42.  NFLAG(I) 

load  code  for  live  loads. 

43.  AL0AD(I,J) 

the  four  loads,  (J  =  1,4)  which  shall  be  the 
target  live  load,  subdivided  into  NJUMPX  incre¬ 
ments  . 

44.  NINC(I) 

NINC(I)  =  0  AL0AD(I,J)  is  the  target  load  to  be 
reached 

NINC(I)  =  1  add  this  load  to  the  previous  target 
load  to  get  the  new  target  load. 

45.  NJUMPX 

the  number  of  increments  into  which  the  next 
load  step  is  to  be  subdivided. 

46.  FC 

ultimate  uniaxial  compressive  strength  of  con¬ 
crete,  entered  as  a  positive  number. 

47.  FTC 

ultimate  uniaxial  tensile  strength  of  concrete, 
entered  as  a  positive  number. 

48.  POISS 

initial  Poisson’s  ratio. 

49.  CEMOD 

initial  uniaxial  elastic  modulus. 

50.  NTYPE 

concrete  constitutive  relation 

NTYPE  =  1  -  linear  uniaxial  concrete  model 

NTYPE  =  3  -  nonlinear  biaxial  concrete  model. 

. 
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51.  NSUB 


the  number  of  subincrements  within  each  load 
increment;  if  NSUB  =  -25,  then  the  program  will 
automatically  choose  the  correct  number  of  sub¬ 
increments  . 


52.  EIO 


the  strain  at  the  maximum  stress  in  uniaxial 
compression,  a  negative  number.  If  a  positive 
value  is  specified,  the  program  will  calculate 
a  value  for  this  variable.  However,  to  use 
this  and  the  following  default,  the  units  of 
input  quantities  must  be  pounds  or  inches. 


53.  EIOT 


the  strain  at  maximum  uniaxial  tension,  a  positive 
number.  If  a  negative  number  is  specified,  the 
program  will  calculate  a  value  for  this  variable. 


A. 2. 3  Load  Codes 

This  section  deals  with  specifying  the  type  of  loads  to  be  applied, 
either  as  forces  or  strains.  Load  codes  are  used  to  specify  dead  loads, 
prestress  loads,  thermal  loads  and  live  loads  and  the  specific  variables 
are  NFLAGD,  NFLAGL ,  NTCODE  and  NFLAG  respectively. 

There  are  four  degrees  of  freedom  for  this  element.  These  degrees 
of  freedom  are  shown  in  Fig.  1.1.  When  a  load  is  to  be  applied  to  the 
wall  segment,  either  a  force  or  axial  strain  must  be  specified  for  degree- 
of-freedom  Nos.  1  and  3,  and  either  a  moment  or  curvature  must  be  specified 
for  degree-of-f reedom  Nos.  2  and  4.  There  are  no  restrictions  for  example, 
that  if  DOF  No.  1  has  a  force  specified  then  DOF  No.  3  must  have  a  force 
specified.  Each  DOF  is  independent  as  to  the  type  of  load  (force  or 
strain)  which  is  specified.  Since  each  DOF  can  be  loaded  in  either  of 
two  ways,  there  are  a  total  of  16  combinations  of  forces  and  strains 
in  which  a  load  vector  can  be  specified. 

To  identify  the  load  combination,  that  is  desired,  a  load  code  is 
used.  The  load  code  is  a  four  digit  number  of  ones  and  zeros  with  each 
digit  corresponding  in  order  to: 


1. 


the  normal  force  or  strain  on  face  1 


2.  the  out-of-plane  moment  or  curvature  on  face  1 

3.  the  normal  force  or  strain  on  face  2 

4.  the  out-of-plane  moment  or  curvature  on  face  2 

In  this  load  code,  a  0  designates  a  strain  or  curvature  and  a  1  designate 
a  force  or  moment. 

For  example,  if  it  was  desired  to  subject  the  wall  segment  to  an 
axial  force  and  moment  on  face  1  and  an  axial  force  and  moment  on  face  2, 
the  load  vector  could  be  described  as 

<  Ni  M]_  N2  M2  > 

and  the  corresponding  load  code  would  be  1111.  If  the  wall  segment  was 

subjected  to  an  axial  force  and  a  curvature  on  face  1,  and  an  axial 

strain  and  moment  on  face  2,  the  load  vector  would  be 

ms 

<  N  i  ^  £  2  M2  > 

and  the  corresponding  load  code  would  be  1001. 

A. 2. 4  Sign  Convention 

(a)  positive  forces  are  tensile 

(b)  positive  moments  cause  tension  on  the  inside  of  the 
element  (closest  to  centres  of  curvature) 

(c)  distances  to  concrete  and  steel  layers  (Z)  are 
measured  from  the  centre  line  along  the  wall  thickness 
of  the  concrete.  Positive  displacements  are  measured 
from  this  centre  line  away  from  the  centres  of  cur- 
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vature  (towards  the  outside  of  the  wall) . 


A. 3  Output 

Program  output  is  read  to  units  6  and  7.  Unit  6  output  begins 
with  an  echo  check  of  the  input  data.  This  is  followed  by  the  results  of 
each  load  increment.  These  results  are  comprised  of  the  specified  load, 
the  stresses  and  strains  in  each  layer,  a  damage  summary  of  the  three 
layers  and  the  overall  results  of  the  specified  load. 

Output  on  unit  7  is  used  mainly  for  program  debugging.  For  each 
load  increment,  the  specified  load  is  first  printed.  Then  after  the 
successful  computation  of  each  load  subincrement,  the  values  of  the  real 
strains,  stresses,  equivalent  uniaxial  strains,  elastic  moduli,  ultimate 
stresses  and  strains  at  ultimate  stresses  are  printed  for  each  concrete 
layer.  This  is  followed  by  the  number  of  iterations  used  to  find  this 
solution.  If  after  30  iterations  no  solution  is  found,  some  additional 
data  is  output,  including  the  trial  real  strains,  the  error  in  forces 
and  moments,  the  total  forces,  the  trial  real  strain  corrections  and  the 
reduction  factor  applied  to  the  strain  corrections.  If  after  121  itera¬ 
tions  no  solution  is  found,  then  the  program  is  terminated.  After  all 
load  subincrements  are  completed,  the  solution  for  the  load  increment 
is  output. 

All  output  solutions  will  be  in  terms  of  units  consistent  with 


those  input. 


' 


APPENDIX  B 

LISTING  OF  COWSAC 


The  author  and  the  University  of  Alberta  will  not  be  responsible 


for  misuse  of  the  following  program  or  for  errors  in  its  listing. 
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APPENDIX  C 

SAMPLE  INPUT  AND  OUTPUT 


C.l  Sample  Input 

The  following  is  the  input  deck  for  example  segment  UD1  as 
discussed  in  Chapter  4. 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 
1 1 
12 
13 
1  4 

15 

16 
1  7 

END  OF 


PROGRAM  TEST  UD1H  EXAMPLE  WALL  SEGMENT 

24.24. ,  12*,  12. ,  1632. ,1632. , 

4,  29600. ,60.  , 

1,1. 1,-9. 3, 

1,  1.  1,  9.  3, 

2,  1. 27,-7. 9, 

2,  1.69,7.9, 

21 11,- 15. 4,- 302.0, 21 .9, -4 6.7, 

2,  2  960  0.,  255.,  1,2, 0,1, 2.  5, 

1,2.  33,0. ,374.  ,  1  , 

2,  2.  33,  0.,  200.  ,  2, 

1  11  1  ,,  6,-12  58.  8,72.  8, -46. 7, 

1  11  1,  .  6,  -12  58.  8,  72.  8,-16  5.  1  , 

1,26,0,1,1, 

26.,  -3. 17,0.  ,.  65E-5,.6  5E-5, 1010, 

1010,6  50.,  19.7  13  2,-532. ,-  3.9332,1  , 

5.  0  0,  .  4  24,.  1  5,  4120.  ,  3  ,-25,  1.,.  000  126, 

FILE 
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C.2  Sample  Output 

The  following  is  a  partial  listing  of  the  output  for  example 
segment  UD1  as  discussed  in  Chapter  4.  This  output  was  generated 
using  the  input  deck  of  the  previous  section  of  this  appendix  and 
includes  the  data  echo  check,  and  the  solutions  to  the  dead  load, 
prestressing,  creep  effects,  temperature  effects  and  the  first  live 


load . 
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